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SUMMARY 

In situations where the experimenter demands the study 
of behaviour of a few factors at the same time and the way these 
factors may interact a form of experimental design is adopted 
known as factorial experiments. 

Due to practical constraints of running an experiment often 
providing homogeneous conditions for all treatments is impossible and 
hence smaller blocks than a full replication of the treatments are 
used. Use of small blocks cause introduction of certain amounts of 
variability to the responser due to heterogenous conditions for the 
blocks. This by itself gives way to a loss of accuracy of certain 
treatment comparisons in the analysis part. 

Previously, work has been done in use of complete and 
incomplete block designs in conjunction with factorial experiments, 
so that the loss in the accuracy of main effects and low order 
interactions is low. 

In this work generalised cyclic designs are introduced and 
work on their properties and mathematical models for their use in 
conjunction with factorial experiments is reviewed. 

Criterion is introduced for distinguishing designs of given 
size, so that the most suitable forms of allocating treatments to blocks 
can be obtained. This criterion is known as efficiency. 

Further, methods of obtaining good designs and past conjectures 
are discussed. 

Finally, tables of cyclic designs and generalised cyclic designs 
with their efficiency values are given. 

Tables of cyclic designs contain information on full sets of 
cyclic designs, within the given range. 



Tables of generalised cyclic designs of order two 
contains the most efficient designs with respect to main effects 
and first order interactions among a large number of designs 
that were run on the EFFICIENCY program. 
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1. INTRODUCTION 

Factorial designs are a useful class of design by which 
the interests of the experimenter demands that effects of different 
factors should be studied simultanious ly . This can be done by 
varying the treatments and testing the treatments under a series of 
differing conditions. Situations where such designs are useful is 
when effects of a few factors have to be determined over a specified 
range and/or it is required to investigate interactions (interrelations) 
between effects of several factors. Some of these factors may be of 
special interest, while others are merely recognizable nuisances. 
The original philosophy and practice of factorial experiments was 
mainly due to work of R.A. Fisher and his book "The design of 
experiments" published in 1935. 

There are several important advantages involved in adopting 
a factorial design rather than the classical arrangement of varying 
only one factor at a time. First of all, one obtains a much more 
general picture of the effect of each factor, because of the range 
of conditions provided by variations in the other factors. Secondly, 
the variety of factor combinations within the experiment facilitates 
the prediction of results likely to be achieved under a particular 
combination of circumstances. Thirdly, if the factors included in 
the experiment have statistically independent effects then as much 
information is available about each factor as would be obtained if the 
whole experiment dealt with only one factor, the rest remaining constant. 
Fourthly, as often happens, the factors are not independent, but have 
effects that are to a greater or lesser degree correlated. In this 
situation only a factorial experiment can provide information about 
the nature of these interactions. This property of the factorial 
experiments will be used to a great extent later on in this project. 



1. 1 INTERACTIONS AND MAIN EFFECTS 

One concept introduced above was the introduction of the 
term interaction, we start by introducing the mathematical model. 
Suppose we are trying to examine the effects of just two factors 
A and B, then we conjecture that observed effect on some experimental 
unit of the ith level of A and jth level of B can be written in the 
form 

y. . = u + a. + b. + c. . + e. .. (1) 
ijk l j ij ljk 

y. ., , is the observed measurement 
ljk' 

u , is an overall average 

a.,b. are relative contributions of the two factors 
i J 

c. . is the interaction contribution 

and e. is the random deviation imposed, 

ljk 

In the analysis part we study the significance of a,b,c 
terms and results are represented in the well known analysis of 
variance tables. 

The affect of a significant interaction can be represented 
in the following figures, for i=l,2 , j=l,2,3 . 




The difference between the two figures is that in the 
case of No Interaction the observed measurements y. . give 
parallel lines, while in Figure 2 the lines are not parallel. 
That is the difference between each level of B does not change in 
Figure 1, no matter what level of A is involved. In Figure 2 
lack of parallelism indicates that apart from main effects some 
additional effects are taking place. This additional effect is 
known as the interaction between the two main effects. 

The above models were only for two main effects present. 
When it is known that several relevant factors are present, the 
case for some kind of factorial design is overwhelming. However 
In situations where a large number of factors are involved the 
magnitude and complexity of the task of design becomes a major problem, 
and often makes a homogenous treatment of a full replication of the 
experiment impossible. Hence we are forced into making a decision 
between involving a large number of treatments or reducing the "efficiency" 
of the design. 

From the elementary form of the model (1) it is possible to 
develop models involving large numbers of factors, block, interactions 
and other sophisticated plans required by practical constraints. 

1.2 FACTORIAL IN COMPLETE BLOCK DESIGNS 

Most types of experimental plans are suitable to be used in 
conjunction with factorial experiments. If the number of treatment 
combinations is not large we can have factorial design in complete 
block designs such as randomized blocks and Latin squares. Complete 
blocks are used with their usual advantages, in such a way as to remove 
effects brought into the experiment due to varying conditions of the 
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experiment. With randomized block designs the grouping of 
experimental units into blocks are used to reduce the error of an 
experiment in one direction only. With the Latin square designs 
the grouping of replications in two directions is employed. 

1.3 FACTORIAL IN INCOMPLETE BLOCK DESIGNS 

For large experiments the major source of error in lack of 
homogenous replications, and the common method of overcoming the 
problem is use of a block which is smaller than a complete replication. 
Class of designs involving block sizes less than a complete replication 
are called incomplete block designs. These designs were introduced by 
F. Yates (1936). In comparison to complete block designs described 
earlier, heterogeneity is eliminated to a greater extent, but the 
reduction in size of a block is achieved by sacrificing all or part of 
the information on certain treatments comparisons. 

1.4 CONFOUNDING 

In factorial experiments use of such blocking will cause the 
sacrifice of accuracy on certain treatment effects, in practice for 
most designs high order interactions are assumed negligible and hence 
main effects and first order interaction comparisons are of prime 
importance. This property can be used so that accuracy of high order 
interactions are sacrificed. 

Using a small block and hence sacrificing high order interaction 
is known as confounding. It will allow precise comparison of main effects 
and important interactions to be made without the involvement of block 
effects . 

It is often possible to construct confounded designs in a few 
replicates in a way that different replicates will cause different 



treatment effects to be confounded, this is known as partial 
confounding. Partial confounding allows some information to be 
retained on all effects and block to block variation is eliminated 
at the same time. 

It is obvious that the best designs are always designs that 
allow a full replication of the treatments to be made in one block. 

As will become clear later we are only concerned with 
multireplicate designs that is designs that will have partial 
confounding. 

1.5 BALANCED INCOMPLETE BLOCK DESIGNS 

One special form of incomplete block designs are balanced 
incomplete block designs, these designs ensure that all pairs of 
treatments are compared with the same precision, this can be achieved 
by retaining balance in the design so that every pair of treatments 
occur together in the same number of blocks and hence all treatment 
comparisons are of equal accuracy. The problem is that these designs 
only exist for a certain number of treatments, block size and number 
of replications. In situations where there are no balanced designs 
further sacrifice of accuracy should be made such that certain treatment 
contrasts have lower accuracy than others. 

1.6 PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS 

As an extension to balanced incomplete block designs, these 
are a group of designs known as partially balanced incomplete block 
designs represented by (PBIB/m), by which, rather than every treatment 
occuring the same number of times, there are m associate classes, and 
hence treatment comparisons can be estimated with m degrees of accuracy. 
As an example group divisible designs are a large and important class of 



partially balanced incomplete block designs with 2 associate classes 
(PBIB/2). t treatments may be divided into m groups each 
containing treatments, where t=m m . Any two treatments in 

the same group are first associates and two treatments in different 
groups are second associates. 

1.7 SOME REFERENCES ON FACTORIAL EXPERIMENTS IN (PBIB/2) 

In the class of partially balanced incomplete block designs 
with two associate classes Kramer and Bradley (1957) and Zelan (1958) 
considered group divisible designs and Bradley, Walpole and Kramer (1960) 
considered Latin square designs with two association schemes. Brenna 
and Kramer (1961) considered the use of factorials in rectangular lattice 
designs. Kurkjian and Zelan (1962,1963) introduced a special calculus 
for the analysis factorial experiments and gave a class of designs having 
an orthogonal factorial structure for two factor experiments, and 
John (1973) has considered use of cyclic designs. 

1.8 CYCLIC DESIGNS 

Cyclic designs are a group of designs belonging to a class of 
incomplete block designs. As the name suggest by the cyclic development 
of the initial block(s) it is possible to obtain all other blocks of the 
design. 

Bose and Nair (1939) have suggested by making the right choice 
of initial block and its cyclic development it is possible to construct 
designs with high degree of balance giving balanced incomplete block 
designs and partially balanced incomplete block designs. 

Cyclic designs as a class in their own right were introduced 
for blocks of size 2 by Kempthorne (1953), Zoellner and Kempthorne (1954) 
and McKeon (1960). Designs for block size 2 were more fully studied 



in two papers by David (1963) and (1965). In these papers it 
was shown how non isomorphic designs of a given size can be 
constructed. These results were extended to general block sizes 
by Wolock (1964), David and Wolock (1965), John (1966). 

Altogether up till now all balanced and partially balanced 
designs considered do not have the necessary flexibility required by 
the factorial experiment. Factorial experiments are multif actor 
experiments where the number of factors involved can be high, and 
in one experiment we are concerned with main effects as well as the 
interactions . 

In the case of balanced data, we perform an analysis of the 
orthogonal treatment combinations. 

In the case of non orthogonal block designs it is possible 

to perform an orthogonal analysis of the main effects and interactions 

provided certain conditions are present. These conditions were shown 

by Cotter, John and Smith (1973). It is difficult to explain this 

result at this stage without the necessary mathematics. However, 

it requires that the design can be represented by a matrix A, that has 

a block circulant property. 

That is for an n factor design with levels m, ,m_. . . .m 

1 2 n 



m. 



A = 



°1 



and 



A. = 
1 

for i=l, . . . ,m^ 

and so on for A^^... such that the last level of blocking has 

a square matrix of order m with row and column sums all equal. 

n 

Designs that have this property are said to have a factorial 
structure. Cyclic designs mentioned previously have a property that 
every n cyclic designs of size k can be combined to give a 
generalised cyclic design of order n represented by GC/n . 
John (1973) has sh own that this generalised cyclic design have 
a factorial structure. 

Later on in the thesis when the mathematical theory is more 
developed results of Cotter, John and Smith (1973) and John (1973) will 
be shown. However at this stage these results were given only to 
show the relation between cyclic designs, generalised cyclic designs 
and factorial structures. 



A li 

A . 
. m 2 1 



A 2i 



2i 



A li 
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2.1 GENERALISED CYCLIC DESIGNS (GC/n) 



Generalised cyclic designs belong to the class of 
incomplete block designs. For an n factor generalised design 
we have an n tuple a^.-.a^ set, representing a particular 
treatment . 

a. is an integer from o to m.-l where m. is the 
1 11 

number of levels associated with the factor i (i=l,...,n) . 

EXAMPLE 

We have two factors a^ and a^ each at 3 and 2 levels 
respectively. Then levels of a^ is represented by 0,1 or 2 and 
a^ is represented by or 1. The complete set of all the possible 
treatments of this 3 by 2 design is 

00, 01, 10, 11, 20, 21 

This ordering of treatments in such a way that the right most digit 
in varying fastest is called LEXICOGRAPHICAL ordering. Now for a 
block size of k=3 we can have (00, 10, 21) as one possible block 
having three of the treatments occurring only once. It is clear 
that considering the first digits we have a cyclic design (GC/1) 
for the factor a^ that is (0, 1, 2), and considering second 
digits we have a cyclic design for a^ that is (0, 0, 1). 

By cyclic development of the above 2nd order generalised 
cyclic design we obtain the full set of the design with 6 blocks. 

(00, 10, 21) (01, 11, 20) 

(10, 20, 01) (11, 21, 00) 

(20, 00, 11) (21, 01, 10) 

Cyclic development is done in such a way that on columns we cycle a. 
and on rows we cycle . 
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The total number of treatments is 6 which is given by 

m l Xm 2 " For a S eneral case wi - th n factors the total number of 
treatments in 

t = m, m ... m 

12 n 

We follow the same representation of the parameters of the GC/n design 
as the one used by John (1973). 

(t,k,r; m , . . . ,m ) where t = total number of treatments 
1 n 

m. = number of levels associated with factor a. 
1 l 

n = number of factors 

k = block size 

r = number of replications. 

The 3x2 design with initial block (00, 10, 21) is therefore represented 
by (6, 3, 3; 3, 2). 

It is clear that there may be many initial blocks with their 
own full set of blocks that have the above representation of the 
parameters. However they may not be of equal interest as far as 
comparison of main effect and interactions are concerned. 

All treatments in the above design occur an equal number of 
times, but we note that the number of times two treatments occur together 
in the same block is not the same. 00 and 11 occur together in the 
same block twice but 00 and 01 do not occur together in any of the blocks. 

The manner in which two treatments occur together in blocks 
is if great importance particularly later when we consider methods by 
which one initial block is to be preferred to another. 

2.2 PROPERTIES OF GENERALISED CYCLIC DESIGNS 

(A) Nonisomorphic sets: Other members of a set generated 
by a particular initial block have the same properties as far as accuracy 
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of comparison of main effects and interactions are concerned. In 

such situations what actually is taking place is that by generating 

the other blocks of an initial block we are relabelling the treatments 

differently. Among these blocks only one is sufficient to give all 

the necessary information we require on the design. In these cases 

we obtain the block giving the lowest numerical value among all 

equivalent blocks so that we have a unique representation of the 

designs and hence obtain a group of non-isomorphic sets. 

We define an operator R(h„ ...h ) such that operation on an 

In 

n-tuple treatment a, ...a changes it to an n tuple treatment b n ...b 

In r 1 i 

(a, ...a )R(h ...h ) - b,...b (2) 

1 n 1 n In 

or a,R(hJ . . . a R(h ) = b. . . .b (3) 

11 n n In 

where h. is a non zero number less than m. relatively prime to 
l i 3 r 

m. . R(h„...h ) can take many different values and operation on 
l 1 n 

initial block (a, ...a ) can produce the other blocks of an equivalent 
In 

class. However it must be noted that many of the blocks produced by 

this operation may not even belong to the full set of (a^. . .a^) and 

so they can produce other blocks by a cyclic rotation of the treatments. 

The way in which h^ produces a relabelling of a^ is as follows, 

h. is multiplied by a. and then mod m. is applied when h. * a. is 
l l i ii 

greater than m. . 

Consider a set of sizes (12, 4, 4; 4, 3) as defined above. 
All possible values for RCh^jh^) are 

R(l,2), R(3,l), R(3,2) 

then 

(00,01, 11, 20)R(3,1) = (00,01,31,20) 

Applying cyclic rotations on either of these two equivalent blocks we 
obtain 12, corresponding equivalent blocks. It can be shown that by 



the above R(h ^b^) operations and cyclic developments, all the 
following blocks are equivalent. 

(00,01,11,20), (00,01,20,31), (00,01,10,21), (00,01,21,30) . 

Block (00,01,10,21) has the lowest numerical values and is 
therefore the initial block. 

If the elimination of the equivalent designs are performed 
systematically it is a very efficient method of reducing the total 
number of blocks to a few intial blocks. 

(B) Fractional sets: 

Consider a set of size (16,4,4; 4,4) with the initial block 

(00 11 22 33) 
Performing cyclic rotations on the two factors we obtain 



(00 


11 


22 


33) 


(10 


21 


32) 


03) 


(20 


31 02 


13) 


(30 


01 


12 


23) 


(01 


12 


23 


30) 


(11 


22 


33 


00) 


(21 


32 03 


10) 


(31 


02 


13 


20) 


(02 


13 


20 


31) 


(12 


23 


30 


01) 


(22 


33 00 


11) 


(32 


03 


10 


21) 


(03 


10 


21 


32) 


(13 


20 


31 


02) 


(23 


30 01 


12) 


(33 


00 


11 


22) 






All 


the 


4 blocks 


in 


any 


row or 


co lumn 


are 


identical with 


4 



blocks in any other row or column. Therefore the full set is not 
composed of 16 blocks as may be expected but only 4 unique blocks. 
We refer to this set as a 1/4 fractional set. 



If t and m. have a common devisor p. for each particular 
i l 

i then a fractional set of size (t,k,k/p; m, m ) exists for p 

In 

where p = p p . . .p 
1 / n 

In the above example P=4=2.2 and so the fractional set is 
of size (16,4,1; 4,4). The full set is equivalent to replicating 
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the fractional set 4 times. For every fractional set of size 

(t,k,k/p; m . . .m ), there exists an equivalent full set represented 
in 

m. m 

by (t/p,k/p,k/p; — — ). 

p l p n 

In the example full set is of size (4,1,1;2,2) given by 
(00), (01), (10), (11). 

Generalised cyclic designs by their high degree of balance 
and also by the way they are able to restrict the number of possible 
permutations of treatments in a block (Property A) are attractive 
designs for situations where a number of treatments exceeds the 
maximum sensible value for a block size. They are very flexible 
designs and due to property B it is always possible to form some 
kind of design for given values of t,k,r and m^ . 

Ease of representation is another advantage of these designs. 
The initial block contains all the necessary information for a design 
of given size, the rest of the design can be obtained from the circulant 
property of the design. This is a major advantage for analysis or 
situations where a comparative study of designs are involved. 

Any GC/n design can be easily extended to GC/n+1 designs 
this dimensional property of generalised cyclic designs makes them 
particularly suitable to be used in conjunction with factorial and 
split plot designs, where some interactions and main effects are more 
important than others. 

2.3 LINEAR MODELS FOR GENERALISED CYCLIC DESIGNS 

We have assumed that all treatments within a block are performed 
under similar conditions, this follows that when an observation is made 
on a particular treatment certain amounts of variability will be due to 
block effects. 
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Model for representing the above block variation is as 

follows 

y. . = u + t. + b . + e. . (4) 



where y. . are the observations 

u is the general mean 

t^'s are the treatment effect 

b.'s are the block effects 
J 

and e^.'s are the uncorrelated random deviations with 
mean zero and constant variance S 2 . 



Alternatively (4) can be written in matrix form 

y»(U t y(utb)' +e (5) 

t and b are vectors of elements t. and b. respectively, and 

l l r 1 ' 

(1 X X^) is the design matrix. 

X^_ and X^ are matrices representing presence or absence of a 
particular treatment effect or block effect. 
X fc and X^ have the following properties. 

(1) Summing rows of X fc or X^ give number of times treatments 
occur in design or number of units in each block respectively, 
that is vectors r and K . 

(2) X. X/ is called the incidence matrix N with elements n. . 

t b 

n. . is the number of times treatment i occurs in block i . 

(3) Where y is the vector of observations. 

y = T vector of treatment totals, 
y = B vector of block totals, 
l'y = G Grand total. 
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Now 



(i x t x t x b ) 



I'l 


1X» 
t 


1X b 
a 


x'l 

t 


X X 
t t 




D 


o t 


b b 


n 


r 1 


K* 




d 


N 


r 


r 


K 




K d 



(6) 



Vectors K and r are row sums of X^ and X^ . Subscript d 
represents diagonal representation of a vector. 



(1 X t X b )'(y) = (G T B)' 
In terms of the usual normal equations X'y = X'X we write 



(1 X fc X b )*(y) = (1 



X 



giving 



nu + r't + K'S = G 

ru + r d t + Nd" = T 

KU + u't + Kb = B 

Eliminating block effects b from the equation, by calculating 



(10) - NK d (11) 



we obtain 



(r-NK 'HOu + (r d -NK _d N')t + (N-NK V^b 

since (r-NK _d K) = . 
We have 

(r^NirVKt) = T - NK _d B 



T - NK d B 



Let 



A = (r d -NK , matrix representing design 

,-d_ 



(7) 



X b )*(l X fc X fa ) (8) 



(9) 

(10) 

(11) 



(12) 



(13) 



and Q = T - NK B , vector of adjusted treatment totals 
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giving 

A t = Q (14) 

(14) is the reduced normal equation. We let W = (A+J) 

be a generalised inverse of A where J is a txt matrix of 

l's. 

Then t = WQ (15) 

Due to special structure of the cyclic designs we have equal block 
sizes and equal replications of treatments. 

This makes r^ and diagonal matrices such that constant 

r and K are multiplied by the identity matrix. 

This simplifies A and Q as follows 

A = rl - (1/K) NN' (16) 

Q = T - (1/K) NB (17) 

Furhter for generalised cyclic designs we can relate 
structure of NN' to structure of A . 

EXAMPLE 

Design of size (12,4,4;4,3) with initial block (00,11,21,32) 



is 



(00, 11, 21, 32) (01, 12, 22, 30) (02, 10, 20, 31) 

(10, 21, 31, 02) (11, 22, 32, 00) (12, 20, 30, 01) 

(20, 31, 01, 12) (21, 32, 02, 10) (22, 30, 00, 11) 

(30, 01, 11, 22) (31, 02, 12, 20) (32, 00, 10, 21) 
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NN' 





00 


01 


02 


10 


11 


12 


20 


21 


22 


30 


31 


3? 


00 


~4 








1 


3 








2 


2 


1 




3 


01 





4 








1 


3 


2 





2 


3 


i 





02 








4 


3 





1 


2 


2 








3 


1 


10 


1 





3 


4 


o 


o 


1 


3 







o 


9 


11 


3 


1 








4 








1 


3 


2 





2 


12 





3 


1 








4 


3 





1 


2 


2 





20 





2 


2 


1 





3 


4 








1 


3 





21 


2 





2 


3 


1 








4 








1 


3 


22 


2 


2 








3 


1 








4 


3 





1 


30 


1 


3 








2 


2 


1 





3 


4 








31 





1 


3 


2 





2 


3 


1 








4 





32 


3 





1 


2 


2 








3 


1 








4 



(18) 



Elements of NN' in this particular case show the number 
of times each treatment occurs with other treatments. Due to this 
property NN' matrix is called the concurrence matrix. Also note 
that treatments are ordered in LEXICOGRAPHICAL ORDERING. 

NN' has a block circulant structure which is similar to A , 

(A+J) and (A+J) . Hence by having one row or column of any of the 

above matrices it is possible to generate the whole matrix. 

We define R^ to be a circulant matrix with first row 
j 

having a 1 in column i and zeros elsewhere. 

nu represent size of the matrix. 

C. . is the value in the jth column of the ith 

block of NN' . 

So that for (18) we have 



C ll = 4 



C 23 = ° 



C 41 = 1 
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Then for the example (Recall 18) 



4 3 . . 

NN f = I I C . (R 1 . id) 
i=l j=l 



(19) 



For a 3 dimensional generalised cyclic design we have 



m l m 2 m 3 



NN 



' - I I I C. .. (R 1 R 2 R 3 , 

i>i i^i i:=i ^ k m l m 2 m 3 ) 



(20) 



For GC/n designs A is given by 

m, m 
1 n 

A = I •'• I a ,- i £ 

i =1 i =1 12 3 
1 n 



where 



n 



L i x l L 2 
n 8R J =R J -0R Z 

j=l m i m i m ? 



n 8 R J 
j-1 



m. 



3 J 



. . 8 R 



n 



m 



n 



(21) 



designs , 
namely 



Since GC/n designs are partially balanced incomplete block 
It follows that W has the same pattern as A , (John 1973) 



m, m 



n 



W = 1 ... I 



w. 
l. 



n 



n ® r 

m. 



j-1 



(22) 
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3.1 EXTENSION OF GC/n DESIGNS TO FACTORIAL MODELS 



So far in the considerations of the model building we have 
been concerned with structure of the blocks used. 
Recall (4) 

V. . = u + t. + b . + e. . 

(i=l,2,...t) 
(j=l,...k) 

Now consider the treatment part of the model, that is the 
t^'s . As we mentioned before in factorial experiments we conjecture 
that value of every observed treatment can be partitioned into part 
due to its main effects, its 1st order interactions and so on for 
higher order interactions. Therefore the model can be written for 
i=l. . . t as 

n n 

t. = I a (i ) + I a (i ,i )+ ... + a ) (23) 

i tl P P p <q=l P q 12... n 1 n 

a refers to main effect of factor p and a to interaction 

p pq 

between factors p and q etc. i„...i are levels of each 

In 

factor corresponding to treatment combination i . 
Recall (7) 

(1 x t x b )'(y) = (GIB)' 

T is the vector of treatment totals. In lexicographical order we 
write T . 

t' = ct t t t 

L u OO...O 0...1 0...2 0...(m-l) 

n 

T 00...10 ... T 0...(m -1) (24) 

n 



7~ t\ /■ in ^(m.-l). . . (m -1) 

(m.j-1) (m -1) . . .0 1 n 
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Similarly we represent t as the vector of treatment effects 
f - (t, 



'0..0 0..01 



t (m 1 -l). .. (m -l) 5 
i n 



Note the subscripts of t^ so far referred to decimal numbers 

i=l...t but now they refer to corresponding number between 
1 to t , expressed by levels of the n tuple treatments. 



EXAMPLE 



For n=2 with m^=4 and m2=3 , we write the vector of 



treatment effects. 



t (t 00 t 01 ^2 t 10 t ll t l2 t 2l t 22 t 30 t 31 t 32 ) 



The average effect of first factor at the ith level is 



1 

3 



C i0 + hi + hi 



for i=0,l,2 

To test main effect of first factor, we test 



By putting 







3 


t 


12 L 

iZ i=0 


I 

j=o 


can be 


expressed 


as 






= 


H o ; ' 


J- • • • 


= 




C 2 " t 


= 



= Overall mean 



In matrix notation the null hypothesis is 



c 10 t = 
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where 



10 



12 



2 2 2 2 -1 -1 -1 -1 -1 -1 -l -1 
-1 -1 -1 -1 2 2 2 2 -1 -1 -1 -1 
-1 -1 -1 -1 -1 -1 -1 -1 2 2 2 2 



J_ 
12 



2 
-1 
-1 



-1 -1 
2 -1 
-1 2 



(1111) 



12 < 3I 3 " V 



We define 



— r 1 « f0 
12 C 3 M C 4 



C = < 
m. 
J 



m.I .-J . 



1* . 



if x.=l 
J 

if x.=0 
J 



The subscript in is a binary number and indicates which 

factorial effect is involved. 
For a 3 factor experiment 



ilOO 



^010 



represents main effect of 1st factor 
represents main effect of 2nd factor 



C 110 represents interaction between 1st and 2nd factors 



For a general case of n factor experiment we write the 



contrast matrix C as 



x x x x 
C X = C 1 n = C 1 6 C 2 
m l m 2 



. 8 C 



m 



(25) 
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where 



x . 

C J 
m. 



m.I -J if x.=l 

J m. m. j 
J J 



1' 



m 



if x.=0 
J 



(26) 



A generalised interaction may be represented by a x . 
Kurkjian and Zelan (1962) have shown that the estimates of the 

• • X 

generalised interaction a is 



a x = e x t 

Recall (15) for the GC/n design 

t = WQ 

Now for factorial experiments we have 

a X = C* WQ 
The following results follow 

X 

(1) Assuming a is estimable then 

a ) = C W C S 2 



(27) 



(28) 



(29) 



(2) As mentioned at the end of introduction GC/n designs have 
factorial structures. 

We have factorial structure if the adjusted treatment sum 

of squares from the block part of the model can be partitioned 

orthogonally into main effect and interaction sum of squares. 

Where an orthogonal factorial structure is defined to be 

C X W C 7 = for all x=fy 

Result of Cotter, John and Smith (1973) for presence of factorial 
structures requires 



m„ 



m 



t 



i =1 
1 P 



n e r 
j-i 



m. 



W. 



i n . . .1 

1 P 



(30) 
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where W. ...i is a square matrix of order m with all column 
sums and row sums equal. 

For GC/n designs, recall (22) 



W= J . 
i. 



w. 



i In 
n 



1 . 

n r J 



(22) and (30) have equivalent structures if 



m 



W 



i i " I 

1 P L =1 

n 



t l 

R n 8 w. 
m i , »«i 
n In 



(31) 



Since W. . in a circulant matrix all row sums and column sums 

1 _ • • • 1 

1 P 

are equal and therefore GC/n designs have orthogonal factorial structures, 



(3) Supposing orthogonality of factorial structures 

c* w c y - - ' 



C W x=y 
x+y 



giving C = C W A as the estimability condition. 



(4) Sum of square due to the hypothesis H 



(32) 



where 



H Q : t = is 



S.S(H Q ) = ((f t)*(C? W C X ')~(C X t) 

= (C X T)'(C X t) 

T is the vector of unadjusted treatment totals, 
t is the vector of treatment effects. 



(33) 



3.2 BALANCE 

If a design has balance with respect to all main effects 
and interactions it is said to have factorial-balance. The analysis 
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of designs with factorial balance is straight forward. However 
it is not always possible to construct designs with factorial 
balance . 

We let h to be an eigenvector of W . Then h't is 
called a BASIC CONTRAST. The importance of basic contrasts is based 
on the independence of the estimators of their treatment effects. 
Suppose there are two basic contrasts b| t and h^t then if they 
are orthogonal that is h^h 2 =0 , the covariance between estimators 
h| t and h^ t is zero and hence they are independent. 

For a generalised interaction h can be represented by 

n x i 

h = H 8 h. J (34) 



j-1 



x . 

where h.-' is a m.xl vector 



and h.^ = • 



x 

1.* 
J 



contrast vector if x.=l 

J 

vector of l's if x.=0 

J 



Now using (33) we obtain the factorial sum of square for the 
basic contrasts. 



H Q : h' t - 

SS(H Q ) = (h' t)' (h W h')" (h' t) 

= (h' t)' (h» t)/L(h' h) (35) 



L(h' h) is a scalar 
But h' t = h 1 W T From (13) and (14) 

(Since t = W Q , with the correction factor zero) 

and h' W T = L h' T 

where L is an eigenvalue of W 
and h is an eigenvector of W . 
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Substituting in (35) we have 

•) - (L hj T)' (L h' T) 

SS( V lTI) 

L(h' T) ' (h* T) 
" (hMT) (36) 

If none of the eigenvalues of corresponding eigenvectors 
of W are zero. Then it would be possible to find an orthogonal 
decomposition of the factorial sum of squares into single degrees of 
freedom component (34, 36) . Given that one or more of the eigenvalues 
are zero then it shows that the corresponding main effect or interaction 
is not estimable. However we have already assumed that all main 
effects and interactions are estimable. 

A main affect is balanced if all its corresponding treatment 
differences can be estimated with the same accuracy. 

Now an interaction is balanced if all eigenvalues corresponding 
to subsets of h are equal. This implies that all linear combinations 
of the eigenvectors have the same eigenvalue. If all eigenvalues 
corresponding to a set of eigenvectors spanning the contrast space (for 
a generalised interaction) are distinct, then any linear combination 
of these eigenvectors will give a unique set with a distinct eigenvalue. 

Analysis for a design with factorial balance is quite simple. 
The sum of squares can be obtained from SS(H ) = L(h' T) 2 /h' h 
and the variance from S 2 L(h' h) . 

For an S factor interaction we let 

h = h* 1 ( S £n) (37) 

(without loss of generality we consider the first s factors) . 

And 1 is a vector of size m , x... x m 

s+1 n 
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A X X 

Then h = h 1 ... 8 h s (x, to x = 1) 

1 s Is 

X^ X 

and h = h. 8 ... 8h n (x , to x = 0) 

1 n s + 1 n ' 

John 1979, b has shown that if h is an eigenvector of 

NN' with eigenvalue L , then h is an eigenvector of N N' 

s s 

with the same eigenvalue. N g is the incidence matrix of a GC/s 
design and the initial block is obtained by taking the 1st s 
elements of the treatment combinations in the initial block of the 
GC/n design. 

EXAMPLE For n=3 , s=2 . 

An intial block for GC/3 is 000 Oil 101 210 
for GC/2 is 00 01 10 21 . 

It is clear that for any GC/n design N g can be obtained 
simply by summing certain elements of N . In this case where the 
interaction of the 1st s factors is considered it is the collapsing 
of the design matrix from the (s+1) th to n th blocking levels. 

For a GC/n design, if a GC/s design has factorial balance 
then all lower order interactions involving these s factors have 
also balance. Hence good designs can be built up from balanced GC/1 
and GC/2 designs. 

3.3 EFFICIENCY 

When we are dealing with a comparative study of experimental 
designs, for a given size there may be a large number of designs to 
choose from. Therefore we need a criterion to distinguish between 
the designs and choose a design that suits our demands. A common 
criterion is known as EFFICIENCY. It is a method of comparing the 
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variance of possible designs with one standard design. This 
standard design is taken to be a randomised block design. 



EFFICIENCY = 



Average variance of estimated treatment contrasts in 

randomised block . 

Average variance of estimated treatment contrasts in 

a particular design. 



It is assumed that error variances s 2 is the same for 
all designs. The Efficiency factor has always a value less than or 
equal to 1. The reason is that the randomized block design by involving 
a full replication of treatments an equal number of times produces a 
lower bound on the variance. 



EXAMPLE 



For a pairwise comparison of t^ and t^ we have 

var(t^-t,,) = 2 s 2 /r , for a randomised block design. 

For some other design we have 

var(t r t 2 ) = (w n + w 22 - 2w 12 )s 2 
- c l2 s 2 

where w. . are values of matrix W with n=l . 

Then the efficiency factor for comparing 1st and 2nd treatments 

is E 12 

„ _ 2 s 2 /r 2 



'12 c 12 s' z rc i2 



giving the overal efficiency factor 

E = ^ 
c 

where 



2 s 2 H c /t(t-l) 
i>j J 
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Efficiency for a factorially balanced design . 

For a randomised block design we have W = l/r(I) 

For a factorially balanced design variance is L(h'h)s 2 . 



Then using (35) we obtain that efficiency for a factorially balanced 
design is (rL) ^ . 

Now in a GC/n design with n factors John (1973) has shown 
the efficiency of the main effect of the ith factor in 



m.-l 

EGO --=: 



rc 
m. 

c = I" l\ 
h=2 n 



(38) 
(39) 



m. 
i 



L = r - (1/k) I n. cos((h-l)(A-l) — ) (h=2...m.) 

11 £=1 1 m i 



(40) 



where n^ are the elements of the 1st row of NN' matrix. 



In a generalised interaction we have (recall 29), 



var(C X t) = C X W C X 



W = (l/r)I is a randomised block design and hence 



E(k) ■ traC£ < C * gX 

r trace (C x W C X 



(41) 



Due to orthogonality of treatment contrasts C W C 1 



= C X W 



(Recall 32) giving 



E(k) 



= trace (C X ) 
r trace (C X W) 

n x - 
.n. (m.-l) J 
= izi J_ 

r trace (C X W) 



(42) 
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John (19 79, a) has shown that, overall efficiency in 
GC/n designs can be related to the efficiency of lower order 
interaction.- The average efficiency for generalised interacti 
C X t is 



on 



E(x) = V / ■ 
x 



V E 1 - J V 

X S L 



x(s) 



x(s) 



E 1 [x(s)] 



(43) 



where E is the overall efficiency for a GC/s desi 



E[k(s)] is the efficiency of main effects up to s factor 
interaction depending on the binary value of (s) 
\ = C™! -1 ) (n^-l) • • • > degree of freedom associated 



with C t . 



V k(s) in de § rees of freedom associated with generalised 



interaction t . 



Using the above result it can be shown that for n=l 
since there are no lower order efficiency factors 

,-1 



I V / E [x(s)] = 
x(s) X(s} 



givxng 



E(x) = 



(m-1) 



(m-l)E 



-1 



Hence overall efficiency in equivalent to the main effect efficiency. 

The importance of this result is however based on the way 
higher order interactions can be evaluated. 

Having the main effect efficiencies for an n factor design 
we can obtain the two factor interactions as follows 

E(110...0) = (m 1 -l)(m 2 -l)/[( I n i m 2 -l)E- 1 - ^-DE^ClO. .OMtyDE^COlQ, . .0)] 

(44) 
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This is a very convenient result for factorial experiments. 
Main affects have highest priority for efficiency and so on for 
lower order interactions. This result shows that using GC/n 
designs it is possible to concentrate the effort of design making 
on obtaining n sets of cyclic designs for the estimation of the 
main effects and then proceed to higher order interactions. This 
result is clearly similar to the final result obtained for balance. 

EXAMPLE 

Consider a design of size (12,4,4; 3,2,2) with the initial 

block 

(000 Oil 101 210) 
For the main effects of the above design we have 

E(100) = 0.94 E(010) = 1.00 E(001) = 1.00 

The above values for the efficiency of the main effects are satisfactory 
and now we consider high order interactions 

E(011) = 1.00 E(110) = 0.81 E(101) = 0.81 E(lll) = 0.44 

We obtain high values for 2 factor interactions. The value 
obtained for the 2nd order interaction is low, however it may be of little 
interest considering (43) it is clear that possibly we can increase the 
efficiency of the 2nd order interaction at the expense of the efficiency 
of 1st order interactions. 
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4.1 CONSTRUCTION OF GOOD DESIGNS 

We have two criteria for distinguishing designs, balance 
and efficiency. The final results of both efficiency and balance 
are similar in that with both criteria, designs have to be built 
for good main effects and then designs with good main effect efficiency 
or balance can be combined in such a way as to give high efficiency and 
balance for the 2 factor interactions. 

From now on the purpose of this work is to construct a set 
of cyclic designs for different values of k and m such that a table 
of relatively efficient designs can be obtained. Once a set of efficient 
cyclic designs is available, for different values of k every two cyclic 
designs will be combined to obtain efficient GC/2 designs for different 
values of m^ and m^ . 

4.2 CONSTRUCTION OF GC/1 DESIGNS 

As was defined previously, size of a GC/n design is represented 
by (t,k,r jm. , . » . »m ) . Number of possible permutations of all possible 
blocks can be obtained from the above parameters. 

Primarily we are concerned with the case of n=l . A block 
size of size k and number of treatments of m gives m k possible blocks. 
However since the position of treatments in a block is not important, we 
consider only the blocks that have their treatments in ascending order. 
This reduces the total number of blocks to (m)(k-l)'. 

It is clear that number of blocks obtained in a set of m(k-l)'. 
blocks is very large. Therefore we use a futher method to reduce the 
number of blocks. Expanding the efficiency formula for the main effects 
(40) 

1^ = r - (1/k) I n £ cos[(h-l)a-l)e] 
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For h=2 writing values of (h-1) (£-1) we have 

1 2 3 .... m.-l 

l 

For h=3 we have o 2 4 6 2(m -1) 

: 1 

For h^ we have (nu-1) 2(m.-l) (m.-l) 2 

if the above values are greater than m. we have to take mode m 

1 i 

so that it is representing one of the treatments. Ignoring the first 
column of zeros the rest of the above values can be represented by a 
randomized block allocation of the m^-1 treatments. 



EXAMPLE For m=5 , (n-1) (£-1) values after taking mod m. can be 
represented by 

h=2 1 2 3 4 

h=3 2 4 1 3 

h=4 3 1 4 2 

h=5 4 3 2 1 

Our intention is to maximize efficiency for the main effect. 

i.e. Minimize C 

and Maximize each 

and Minimize n^cos [(n-1) (£-1) e] 

n 1 is always the greatest contributor to values of . n is 

multiplied by the value of cos(0)=l. While n ...n are multiples 

^ m r 

of values less than 1. 

For a cyclic representation of NN' , the n^^ values are 
in fact the diagonal elements of this matrix, and in binary designs 
(see later) they represent the number of times each treatment is 
replicated. 
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For a general case 



m-1 

j=0 



(45) 



where f(j) are values of the first column of N 

note j+£-l is reduced by mod where necessary. 

and N is a matrix that represents the number of times 
each treatment occurs in each block. 



EXAMPLE 



design 



(00123) 
(11234) 
(22340) 
(33401) 
(44012) 



k=5, m=5 and initial block 00123 produces the following 



The matrix N is 



Block 

1 
2 
3 
4 
5 



Treatments 12 3 4 



2 

1 
1 
1 



1 
2 

1 
1 



1 
1 

2 

1 



1 
1 
1 
2 





1 
1 
1 
2 



is equal to multiplication of any column of N by its own transpose. 
It is minimized if f(i) are equal or as equal as possible. This 
criterion is that the best designs are the ones that give an as equal 
as possible replications of treatments within a block, it reduces the 
total number of blocks to a reasonable amount that can be written in 
a simple way. This set obtained however, contains many equivalent 
designs and theory regarding non-isomorphic sets can be used to reduce 
these designs further, to a few initial blocks. 

Designs having an N matrix with values zero and one are 
called binary designs. One major difference between binary and non-binary 
designs as far as we are concerned is that in binary designs NN' matrix 
(concurrence matrix) represents the number of times two treatments occur 



- 34 - 



together (Recall 18) . The values of the concurrence matrix are 
very important for the efficiency of the designs, and it has been 
conjectured previously by Williams and Patterson (1978), that the 
most efficient designs are to be found in the class of designs where 
the off diagonal elements of the concurrence matrix do not differ 
by more than one. 

EXAMPLE 

Consider design of size k=4, m=5 . 

BINARY: A binary design for the above size has the initial 
block (0123), giving the first column of N , 
(1,1,1,1,0). Then the first colum of NN' using 
(45) is (4,3,3,3,3). It is clear that if a full 
replication of the design is obtained, off diagonal 
elements of (4,3,3,3,3) represent the number of times 
appears with the other treatments. 

NON BINARY: For the same size a non binary design has initial block 
as (0012), giving the first column of N , (2,1,1,0,0) 
and first column of NN' as (6,3,2,2,3) 

The conjecture made by Williams and Patterson (1978) is 
effective for binary designs. When we are concerned with non binary 
designs there are many situations where the off diagonal elements do 
not differ by more than one. However the diagonal elements of the 
design can be decisive in choosing the best design. 

Consider a design of size k=4, m=3 and initial blocks 
(0 12) and (0001). The two designs have (6,3,3) and 
(10,3,3) as the first columns of the NN' matrix respectively. 

On considering the expansion of (40) it was shown that the 
diagonal value of the NN' matrix is very important when choosing a design. 
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Now apart from conjecture of Williams and Patterson we require 
that with non binary designs the best designs are the ones that 
give the lowest diagonal value in the NN' matrix. That is 
between (0 12) and (0001); (0012) gives a higher average 
efficiency. 

m 

Binary designs have £ f (i) = K 

i=l 

giving 

rK = n 1 + n. + . . . + n (46) 
12m ' 

However the diagonal element is equal to the number of replications 
(r) and hence 

r(K-l) = n„ + . .. + n (47) 
l m 

(47) is a useful formula to obtain an optimum design with the 
Williams and Patterson conjecture. 

So far in the consideration of efficiency we have dealt 
with the structure of NN' . This matrix is a very good guide to come 
close to an optimum efficiency design. However designs with similar 
values in the structure of the NN' matrix can give different efficiency 
values, and it is necessary to obtain the efficiency values. Calculation 
of the efficiency values for GC/n designs is relatively complex in that 
information on the structure of these designs is used to speed up the 
algorithm. For our calculations we use the program EFFICIENCY, 
(see appendix III) . 

Input to the program is the size of the design and the initial 
block. The output is the average efficiency for the factorial effects 
and the overall efficiency. 

W is a generalised inverse of A (Recall 16) 

A = R I - (1/k) NN' 
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In the program we let 

W = [r I - (1/k) NN' + j]" 1 (48) 

The efficiencys obtained for the efficiency of the 
best cyclic designs are given in the next chapter. Set of all the 
non isomorphic designs are given with their first column of NN' matrix. 
For the best designs the actual efficiency values are also given. 

One quick method for obtaining the first column of the NN' 
matrix (in binary designs) is to produce only the concurrence of the 
zero treatment . 

EXAMPLE 

K=3, m=7, initial block (0 1 2) then 3 blocks with 
treatment zero are (0 1 2) 

(6 1) 
(5 6 0) 

Giving the first column of the NN' matrix as (3,2,1,0,0,1,2) . 
The above blocks are also given in the tables of the cyclic design. 

4.3 CONSTRUCTION OF GC/2 DESIGNS 

Once the best designs for a given value of K are known. 
Different cyclic designs of size K can be combined for values of m^ 
and m, . Depending on the way two cyclic designs are combined the 
information regarding the interaction between the two factors can be 
used. We can have many GC/2 designs when two cyclic designs are 
combined. However as was explained in the last chapter the value of 
the overall efficiency of the design is directly proportional to the 
interaction efficiency. 

John and Dean (1975) and De an and John (19 75) have given 

general mthods for constructing single replicate factorial experiments 

in generalised cyclic designs for both symmetrical and asymmetrical 
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arrangements. This method is based on obtaining one or more 
generators, other treatments within the initial block are obtained 
by multiplying the generator by integers less than m , and taking 
mods where necessary. 

EXAMPLE 

For a 4 2 experiment a 12 generator forms the inital block. 

00, 12, 20, 32 

where (12)x22 = 24 giving 20 

(12)x33 = 36 giving 32 

In the case of multirepli cate designs, the treatments produced 
from the generator are in fact the treatments that have to be avoided 
within the same block. Considering the above initial block, we write 
other blocks with treatment 00 . 



(00 


12 


20 


32) 


(32 


00 


12 


20) 


(20 


32 


00 


12) 


(12 


20 


32 


00) 



which is clearly not a good design. 

Even if we had a better generator for the main effects, say 11, 
the initial block and other blocks with 00 treatment are 



(00 


11 


22 


33) 


(33 


00 


11 


22) 


(22 


33 


00 


11) 


(11 


22 


33 


00) 



and again as far as second order effect is concerned, the values of the 
NN' matrix are far from the optimal, and hence the overall efficiency 
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is low. The first column of the NN' matrix is 

(4,0,0,0,0,4,0,0,0,0,4,0,0,0,0,4) , we have 4, fours for position 
of 00,11,22,33 and all the other values are zero. The sum of the 
values of the first column of NN' is 16. A better design is the one 
that spreads this value more evenly. 

A better design has initial block (00 11 23 32) giving 



00 


11 


23 


32 


33 


00 


12 


21 


21 


32 


00 


13 


12 


23 


31 


00 



Except for 00 none of the treatments within the initial block 
can be generated from the other treatments. The first column of the 
NN' matrix is (4,0,0,0,0,1,2,1,0,2,0,2,0,1,2,1) . 

In fact the above NN' matrix is very good. Since for the 
two main effects it has full average efficiency, that is E (01) =E ( 10) =1 
and for the interaction E(ll)=0.6279 . Giving the overall efficiency 
of 0.7377. 

When considering design of good GC/2 designs we can assume that 
already we have decided on the value of the maximum average efficiency 
of the main effects. And so the problem is not only finding the highest 
overall efficiency for different combinations of two cyclic designs, but 
also considering other cyclic designs with the same maximum efficiency 
that can give a relatively high overall efficiency for the GC/2 design. 

Two GC/1 designs of size K can yield as many as K 2 different 
combinations of GC/2 designs. Further for each cyclic design there 
may also be a few more designs with the same efficiency and hence the 
total number of designs to consider is very large. 

To overcome the problem of reducing the total number of possible 
designs we make the following observation on the structure of the NN' matrix. 
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Consider design of size (30,4,4;6,5) . 

Cyclic design giving the best efficiency is (0123), for both main effects. 
We write all the treatments for the effects serially. 

Effect A 012345012345012 ... 

Effect B 012340123401234 ... 

The worst possible way of combining these designs is to take 





1 


1 
1 


2 
1 


3 
1 


4 


5 1 





1 


2 


3 


4 


12 


giving 
















00 




11 


22 33 






S4 




00 


11 22 






43 




54 


00 11 






32 




43 


54 00 


A better design 


can 


have 


the combinations 





1 


2 




3 


4 5 





1 


2'"' 




'3 


4 


giving 




00 




11 


23 32 






54 




00 


12 21 






42 




53 


00 14 






33 




44 


51 00 



Therefore using the information on the slope of the lines can 
lead us to a good NN' matrix. 

The above method is a relatively easy method of obtaining almost 
the best designs particularly for the high values of K . When using 
this method one should be aware of the point that the method although 
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often leads to very good designs it is dependent only on the value 
of K . 

Consider designs of the following sizes. 



(24,4, 4;7, 4) 



1 2 3^4 5 6 



1 2' 3 1 ^2 3 



(32,4,4;8,4) 



12 3^ 45 6 7 




(36,4,4;9,4) 



012 3^ 456 7801 



T 3 



In fact for many designs of size K=4 it was observed that 

o o °. o 

combining the treatments in the form I j Y lead to best 

designs while for K=5 the best designs are 





The method arrives at good designs once the two cyclic designs 
are known. However it does not provide information for situations where 
the cyclic design is not decided but the maximum efficiency for the main 
effects is known. In such cases we obtain a set of good designs. 

In the beginning of this section we mentioned the difference 
for the initial blocks of single replicate and multireplicate designs. 
That is for multireplicate designs we require that as few as possible 
of the treatments within an initial block can be generated from each 
other. This is a good method of reducing the number of designs if we 
have a set of good designs. 
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The method is based on obtaining good NN' matrices that 
have their elements almost or as equal as possible. However there 
is no reason to suppose that two designs with equivalent forms of 
NN 1 matrix should give the same efficiency values. In fact for 
K=3 differences of order 15% in the efficiency was noted. 

Consider designs of size (54,3,3;9,6) and initial blocks 
(00 11 34), (00 14 31) 



(00 


11 


34) 


(00 


14 


31) 


(85 


00 


23) 


(82 


00 


23) 


(62 


73 


00) 


(65 


73 


00) 



The off diagonal elements of NN' matrix of the above designs 
are either or 1 but the overall efficiency of the first is 0.3128 and 
for the second is 0.4853. 

However the effect of this property is reduced by a great amount 
such that the difference reduces to order of 5% in the case of K=4 
(on the designs that were studied). 

So far all the designs we considered were composed of two binary 
cyclic designs. Now we consider the method for a non binary design. 

Size (12,4,4;4,3) 01230123... 




giving the initial block (00 11 20 32) and overall efficiency of 
(0.7520). 

It was shown that designs of form | ^ ^ J • • • have low 

efficiency for all values of K . In fact the intention is to obtain 

patterns that will give high efficiency and then use the pattern to 
obtain good designs for other sizes. 
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Judging a design purely on the basis of obtaining different 
slopes, may not always give good designs as is shown with the following 
case . 




Now considering the NN' matrix we have 



00 


11 


23 


35 


47 


88 


00 


12 


24 


36 


76 


87 


00 


12 


24 


64 


75 


87 


00 


12 


52 


63 


75 


87 


00 



In fact it is observed that for high values K not only we should choose 
designs with lowest number of similar slopes but also rate of change 
should also be considered. Chapter 5 contains tables of some GC/2 
designs with relatively high efficiencys. 
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5.1 Table of cyclic designs . 

This table contains information on cyclic designs and 
their efficiency within the following range. 

Block size k, 3 to 8. 

Number of treatments m, 2 to 9 . 

The first column of this table gives the values of k and 
number of replications r . As a standard r is taken to be equal 
to k so that the most efficient designs for each block size can be 
obtained. 

The second column gives the number of treatments. 

The first two columns refer to the practical constraints of 
the experiment and these values have to be known. 

For situations where there are more than one initial block 
for value of k and m the 3rd column gives a reference letter to 
distinguish designs. 

By considering the final column of the table for each value 
of k and m the best initial block cna be obtained. However, the 
other initial blocks are given for completeness. 

The fifth column gives the first column of the NN' matrix. 
This column is of theoretical interest. 
Example 

Block size k=4, number of treatments = 9. 
The table gives 4 initial blocks. 

Designs A and D give the lowest values for the average 
efficiency and in fact by considering the NN' matrix (5th column) it is 
clear that these two initial blocks do not provide an optimal allocation 
of the treatments to the blocks. 

Designs B and C, give NN' matrices that have off-diagonal 
elements differing by 1 only and hence give the highest efficiency. 
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The most efficient design is C with initial block 
(0134) and efficiency of 0.8340. 

In the above tables it is possible to obtain initial 
blocks for particular values of m and k in such a way that the 
diagonal elements of the NN' matrix are higher than the ones given 
in the tables, but as was discussed this will reduce the efficiency 
by a great amount. 
Example 

k=4, m=9 

One possible initial block that is not given in the tables 

is (0012). 

This design gives the value of 6 for the diagonal element 
of the NN' matrix and hence has a much lower efficiency. 

Further, in Appendix I for each initial block all other 
blocks containing treatment zero are also given. 
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TABLE 1 



Table of all initial blocks of cyclic designs 
For Block sizes 3 to 8, number of treatments 2 to 9 





m 


Ref 


Initial block 


1st column of NN' 


Efficiency 


k=3 
r=3 


2 
3 
4 
5 
6 

7 
8 

9 


A 
B 
C 
A 
B 
A 
B 
C 
D 
A 
B 
C 


001 
012 
012 
012 
012 
013 
024 
012 
013 
012 
013 
014 
024 
012 
013 
036 


(5,4) 

(3,3,3) 

(3,2,2,2) 

(3,2,1,1,2) 

(3,2,1,0,1,2) 

(3,1,1,2,1,1) 

FRACTIONAL 

(3,2,1,0,0,1,2) 

(3,2,1,0,0,1,2) 

(3,2,1,0,0,0,1,2) 

(3,1,1,1,0,1,1,1) 

(3,1,0,1,2,1,0,1) 

(3,0,2,0,2,0,2) 

(3,2,1,0,0,0,0,1,2) 

(3,1,1,1,0,0,1,1,1) 

FRACTIONAL 


0.8889 

1.000 

0.8889 

0.8148 

0.7435 

0.7843 

0.6833 
0.7778 
0.6312 
0.7467 
0.6914 

0.5862 
0.7224 


k=4 
r=4 


2 
3 
4 
5 
6 

7 

8 

9 


A 
B 
C 
A 
B 
A 
B 
C 
D 
E 
F 
A 
B 
C 
D 


0011 
0012 
0123 
0123 
0123 
0124 
0134 
0123 
0124 
0123 
0124 
0125 
0134 
0145 
0246 
0123 
0124 
0134 
0136 


(8,8) 

(6,5,5) 

(4,4,4,4) 

(4,3,3,3,3) 

(4,3,2,2,2,3) 

(4,2,3,2,3,2) 

FRACTIONAL 

(4,3,2,1,1,2,3) 

(4,2,2,2,2,2,2) 

(4,3,2,1,0,1,2,3) 

(4,2,2,1,2,1,2,2) 

(4,2,1,2,2,2,1,2) 

(4,2,1,2,2,2,1,2) 

FRACTIONAL 

FRACTIONAL 

(4,3,2,1,0,0,1,2,3) 

(4,2,2,1,1,1,1,2,2) 

(4,2,1,2,1,1,2,1,2) 

(4,1,1,3,1,1,3,1,1) 


1.000 

0.9375 

1.000 

0.9375 

0.8937 

0.8937 

0.8517 
0.8750 
0. 8095 
0.8498 
0.8498 
0.8498 

0.7711 
0.8324 
0.8340 
0.8036 




2 
3 
4 
5 
6 
7 




00011 
00112 
00123 
01234 
01234 
01234 


(13,12) 

(9,8,8) 

(7,6,6,6) 

(5,5,5,5,5) 

(5,4,4,4,4,4) 

(5,4,3/3,3,3,4) 


0.9600 

0.9600 

0.9600 

1.000 

0.9600 

0.9306 


k=5 
r=5 


8 
9 


A 
B 
C 
D 
A 
B 
C 
D 


01234 
01235 
01245 
01246 
01234 
01235 
01236 
01346 


(5,4,3,2,2,2,3,4) 

(5,3,3,3,2,3,3,3) 

(5,3,2,3,4,3,2,3) 

(5,2,4,2,4,2,4,2) 

(5,4,3,2,1,1,2,3,4) 

(5,3,3,2,2,2,2,3,2) 

(5,3,2,3,2,2,3,2,3) 

(5,2,2,4,2,2,4,2,2) 


0.9036 
0.9106 
0.9081 
0.8960 
0.8759 
0.8958 
0.8962 
0.8862 
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m 


ret 


Initial block 


1st column of NN' 


Efficiency 




2 




000111 


(18,18) 


1.000 




3 




001122 


(12,12,12) 


1.000 




4 


A 


001123 


(10,9,8,9) 


0.9623 






B 


001223 


FRACTIONAL 






5 




001234 


(8,7,7,7,7) 


0.9722 


k=6 


6 




012345 


(6 6 6 6 6 61 


1 nnn 
JL # uuu 


r=6 


7 




012345 


(6,5,5,5,5,5,5) 


0.9722 




8 


A 


012345 


(6,5,4,4,4,4,4,5) 


0.9511 






B 


012346 


(6,4,5,4,4,4,5,7) 


0.9510 




9 


A 


012345 


(6,5,4,3,3,3,3,4,5) 


0.9320 






B 


012346 


(6,4,4,3,3,4,4,4,4) 


0.9361 




2 




0000111 


(25,24) 


0.9796 




3 




0001122 


(17,16,16) 


0.9796 




4 




0011223 


(13,12,12,12) 


0.9796 




5 




0011234 


(11 10 9 9 10) 


9689 


k=7 


6 




0012345 


(9,8,8,8,8,8,8) 


0.9796 


r=7 


7 




0123456 


(7,7,7,7,7,7,7) 


1.000 




8 




0123456 


(7,6,6,6,6,6,6,6) 


0.9796 




9 


A 


0123456 


(7,6,5,5,5,5,5,5,6) 


0.9636 






B 


0123467 


(7,5,5,6,5,5,6,5,5) 


0.9635 




2 




00001111 


(32,32) 


1.000 




3 




00011122 


(22,21,21) 


0.9844 




4 




00112233 


(16,16,16,16) 


1.000 




5 




00112234 


(14,13,12,12,13) 


0.9763 


k=8 


6 


A 


00112345 


(12,11,10,10,10,11) 


0.9746 


r=8 




B 


00122345 


(12,10,11,10,11,10) 


0.9746 






C 


00123345 


(12,10,10,12,10,10) 


0.9740 




7 




00123456 


(10,9,9,9,9,9,9) 


0.9844 




8 




01234567 


(8,8,8,8,8,8,8,8) 


1.000 




9 




01234567 


(8,7,7,7,7,7,7,7,7) 


0.9844 
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5.2 Tables of GC/2 designs . 

These tables contain a set of GC/2 designs for block 
size k, ranging from 3 to 8, and the number of treatments (m^ and m 2 ) 
for each of the two factors ranging from 2 to 9. 

For each set of k, m^ and m_, one design is given. This 
is the design that produces the highest overall efficiency and hence the 
highest first order interaction efficiency among a large set of designs. 
Priority was given to the main effect efficiencies and hence the only 
GC/2 designs that were considered were the ones giving the maximum 
possible main effect efficiencies (from Table 1). 

The last 4 columns of these tables give the efficiency 
values for the main effects of the first factor, main effect of the 
2nd factor, interaction effect efficiency and the overall efficiency, 
respectively. 
Example 

For block size 6 and number of treatments 7 and 3, the given 
design is (00 11 20 32 42 51) giving; 

0.9722 and 1.00 as efficiency values for the two main effects, 
0.7854 is the interaction effect efficiency and 0.8529 is the overall 
efficiency of the design. 

In the Appendix II for each design in the table 2, k other 
blocks of the same design that can be produced by a cyclic development 
and contain the treatment 00 are also given. 
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TABLE 2A 

Table of selected GC/2 designs for block size 3 



m l 


m 2 


DESIGN 


E(10) 


E(01) 


E(ll) 


E 


2 


2 


00 


11 


00 


. 8889 


. 8889 


. 8889 


. 8889 


3 


2 


00 


10 


21 


1.000 


. 8889 


.5556 


. 7435 


4 


2 


00 


10 


21 


. 8889 


. 8889 


.5333 


.6914 


5 


2 


00 


10 


21 


. 8148 


. 8889 


.5859 


.5998 


6 


2 


00 


10 


31 


.7843 


.8889 


.4040 


.5535 


7 


Z 


00 


10 


31 


. 7778 


o o on 

. 8889 


. 4866 


.6141 


8 


2 


00 


10 


31 


.7467 


.8889 


.5248 


.6292 


9 


2 


00 


1 

X \J 


31 

«J X 


.7224 


.8889 


.5129 


.6116 


3 


3 
















4 


3 


00 


11 


22 


. 8889 


1.000 


.3439 


.4820 


5 


3 


on 


1 1 

X X 


22 


.8148 


1.000 


.2927 


.4088 


6 


3 


Of) 


1 1 

J_ X 


S9 


. 7843 


1.000 


.4796 


.5817 


7 


3 


00 


1 1 

X X 


19 


. 7778 


1.000 


.5051 


.5975 


8 


3 


00 


1 1 

X X 




.7467 


1.000 


. 4585 


.5488 


9 


3 


00 


11 


32 


.7224 


1.000 


.3895 


.4802 


4 


4 


oo 


1 9 


9 1 


.8889 


.8889 


.5333 


.6349 


5 


4 


on 


1 9 
J-Z. 


9 1 


. 8148 


. 8889 


. 4204 


.5159 


6 


4 


00 


1 9 

X Z_ 


— ' X 


.7843 


.8889 


.4912 


.5709 


7 


4 


00 


11 


32 


.7778 


. 8889 


.3420 


.4238 


8 


4 


00 


11 


32 


.7467 


.8889 


.4460 


.5181 


9 


4 


00 


11 


32 


.7224 


.8889 


.4674 


.5319 


5 


5 


00 


12 


21 


. 8148 


.8148 


.5014 


.5752 


6 


5 


00 


11 


32 


.7843 


.8148 


.4379 


.5092 


7 


5 


00 


11 


32 


.7778 


.8148 


.2991 


.3661 


8 


5 


00 


11 


32 


.7467 


.8148 


.2729 


.3336 


9 


5 


00 


11 


32 


.7224 


.8148 


.3591 


.4187 


6 


6 


00 


11 


34 


.7843 


.7843 


.4724 


.5329 


7 


6 


00 


11 


34 


.7778 


.7843 


.4518 


.5094 


s 


6 


00 


11 


34 


7467 








9 


6 


00 


14 


31 


.7224 


.7843 


.4359 


.4853 


7 


7 


00 


13 


31 


.7778 


.7778 


.4481 


.5012 


8 


7 


00 


11 


33 


.7467 


.7778 


.2157 


.2597 


9 


7 


00 


11 


33 


.7224 


.7778 


.1974 


.2367 


8 


8 


00 


13 


61 


.7467 


.7467 


.4317 


.4763 


9 


8 


00 


13 


61 


.7224 


.7467 


.4240 


.4655 


9 


9 


00 


13 


61 


.7224 


.7224 


.4176 


.4561 
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TABLE 2B 



Table of selcted GC/2 designs for block size 4 



m l 


m 2 


DESIGN 


E(10) 


E(01) 


E(ll) 


E 


2 


2 


00 


01 10 


11 


1,000 


1.000 


1.000 


1.000 


3 


2 


00 


01 10 


21 


.9375 


1.000 


.8125 


.8937 


4 


2 


00 


10 21 


31 


1.000 


1.000 


.6000 


.7778 


5 


2 


00 


11 21 


30 


.9375 


1.000 


.6875 


.8119 


6 


2 


00 


10 21 


31 


.8937 


1.000 


.6224 


.7520 


7 


2 


j oo 


10 21 


41 


.8750 


1.000 


.6742 


.7758 


8 


2 


00 


10 21 


41 


.8498 


1.000 


.6897 


. 7737 


9 


2 


00 


10 21 


41 


.8324 


1.000 


. 7052 


. 7743 


3 


3 




m in 


99 


.9375 


.9375 


.7031 


.8036 


4 


3 


00 


11 20 


32 


1.000 


.9375 


.6319 


.7520 


5 


3 


00 


10 21 


32 


.9375 


.9375 


.5996 


.7091 


6 


3 


00 


10 21 


32 


.8937 


.9375 


.6586 


.7420 


7 


3 


00 


10 21 


42 


. 8750 


.9375 


6049 


6937 


8 


3 


00 


10 31 


42 


. 8498 


.9375 


. 6698 


7355 


9 


3 


00 


10 31 


42 


.8340 


.9375 


.6762 


. 7347 


4 


4 


00 


11 23 


32 


1.000 


1.000 


.6279 


.7377 


5 


4 


00 


11 23 


32 


.9375 


1.000 


.6875 


.7689 


6 


4 


00 


11 23 


32 


. 8937 


1.000 


.6586 


. 7332 


7 


4 


00 


11 23 


42 


.8750 


1.000 


.5967 


.6746 


8 


4 


00 


11 32 


43 


. 8498 


1. 000 


.5666 


.6418 


9 


4 


00 


11 23 


42 


.8324 


1.000 


.6060 


.6703 


5 


5 


00 


11 23 


32 


.9375 


.9375 


. 6036 


6849 


6 


5 


00 


11 23 


32 


. 8937 


.9375 


.6561 


. 7188 


7 


5 


00 


11 23 


42 


.8750 


.9375 


.6690 


.7234 


8 


5 


00 


13 32 


41 


. 8498 


.9375 


.6639 


. 7133 


9 


5 


00 


13 32 


41 


.8340 


.9375 


.6575 


.7037 


6 


6 


00 


11 23 


32 


. 8937 


. 8937 


.6643 


. 7169 


7 


6 


00 


11 23 


42 


.8750 


.8937 


.6487 


.6984 


8 


6 


00 


11 23 


42 


.8498 


.8937 


.6539 


.6978 


Q 





00 


13 32 


41 


. o 


QQ T *7 

. oy J / 


. 040/ 


c on c 

. 6o96 


1 


7 


00 


11 24 


42 


.8750 


.8750 


.6470 


.6921 


8 


7 


00 


11 24 


42 


.8498 


.8750 


.6470 


.6877 


9 


7 


00 


11 34 


42 


.8340 


.8750 


.6394 


.6774 


8 


8 


00 


11 34 


43 


.8498 


.8498 


.6332 


.6712 


9 


8 


00 


11 34 


42 


.8340 


.8498 


.6491 


.6820 


9 


9 


00 


11 34 


43 


.8340 


.8340 


.6438 


.6745 
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TABLE 2C 

Table of selected GC/2 designs for block size 5 



m i 


rn 

m 2 


DESIGN 


EflO") 


E ("01 ") 




lit 


2 


2 


00 


00 


01 


10 


11 


. 9600 


. 9600 


.9600 


.9600 


3 


2 


00 


01 


10 


11 


20 


. 9600 


. 9600 


9600 


9600 


4 


2 


00 


01 


10 


11 


21 


.8471 


.9600 


7385 




5 


2 


00 


11 


20 


30 


41 


1.000 


.9600 


.7179 


.8480 


6 


2 


00 


11 


20 


30 


51 


. 9600 


.9600 


. 7561 


. 8552 


7 


2 


00 


11 


20 


30 


41 


.9306 


.9600 


.7525 


.8407 


8 


2 


00 


11 


20 


30 


51 


.9126 


.9600 


.7582 


.8359 


9 


2 


00 


11 


20 


30 


51 


.8958 


.9600 


.7756 


.8380 


3 


3 


00 


01 


10 


11 


22 


96oo 




. Q4UU 




4 


3 


00 


01 


10 


21 


32 


.9600 


.9600 


.7924 


.8907 


5 


3 


00 


10 


21 


31 


42 


1.000 


. 9600 


6081 

i UUU1 


7277 


6 


3 


00 


11 


21 


32 


42 


.9600 


.9600 


.7200 


.8026 


7 


3 


00 


11 


21 


32 


40 


.9306 


.9600 


.7003 


.7793 


8 


3 


00 


11 


21 


32 


60 


.9126 


.9600 


.7171 


. 7856 


9 


3 


00 


11 


21 


32 


60 


.8962 


.9600 


.7445 


.8000 


4 


4 


00 


02 


11 


23 


30 


9600 


9600 


7322 




5 


4 


00 


11 


20 


32 


43 


1.000 


.9600 


.7317 


.8077 


6 


4 


00 


11 


20 


32 


53 


.9600 


.9600 


.7418 


.8055 


7 


4 


00 


11 


20 


32 


43 


.9600 


.9600 


.7431 


. 7990 


8 


4 


00 


11 


23 


32 


50 


.9126 


.9600 


.7646 


.8102 


9 


4 


no 


1 1 

X A, 


23 


32 


50 


.8958 


.9600 


.7622 


.8038 


5 


5 


00 


11 


23 


34 


42 


1.000 


1.000 


.7040 


.7811 


6 


5 


00 


11 


24 


32 


52 


.9600 


1.000 


.7283 


.7909 


7 


5 


00 


11 


23 


34 


42 


.9306 


1.000 


.7263 


.7817 


8 


5 


00 


12 


21 


34 


53 


.9126 


1.000 


.7398 


.7876 


9 


5 


00 


12 


21 


34 


53 


895 8 


1 ooo 


7513 


7Q2S 


6 


6 


00 


11 


23 


34 


42 


.9600 


.9600 


.7273 


.7814 


7 


6 


00 


11 


22 


34 


43 


.9306 


.9600 


.7272 


.7749 


8 


6 


00 


11 


23 


54 


62 


.9081 


.9600 


.7436 


.7837 


9 


6 


00 


11 


23 


34 


52 


.8958 


.9600 


.7069 


.7494 


7 


7 


00 


12 


21 


34 


53 


.9306 


.9306 


.7510 


.7897 


8 


7 


00 


12 


21 


34 


53 


.9126 


.9306 


.7350 


.7718 


9 


7 


00 


12 


21 


34 


53 


.8958 


.9306 


.7449 


. 7768 


8 


8 


00 


11 


23 


36 


62 


.9126 


.9126 


.7353 


.7684 


9 


8 


00 


11 


23 


35 


52 


.8958 


.9126 


.7150 


.7480 


9 


9 


00 


11 


23 


35 


52 


.8958 


.8958 


.7026 


.7343 
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TABLE 2D 

Table of selected GC/2 designs for block size 6 



m 

m i 


m 
2 


DESIGN 


E(10) 


E(01) 


E(ll) 


E 


2 


2 


00 


00 


01 


10 


11 


11 


1.000 


1.000 


.8889 


.9600 


3 


2 


00 


01 


10 


11 


20 


21 


1.000 


1.000 


1.000 


1.000 


4 


2 


00 


01 


10 


11 


ZD 


Q 1 
J 1 


.9623 


1.000 


.9252 


.9510 


5 


2 


00 


01 


10 


20 


31 


41 


.9722 


1.000 


.8432 


.9129 


6 


2 


00 


10 


21 


30 


41 


51 


1.000 


1.000 


.7843 


.8889 


7 


2 


00 


11 


21 


30 


41 


50 


.9722 


1.000 


.7856 


.8779 


8 


2 


00 


11 


21 


30 


41 


60 


.9510 


1.000 


. 7723 


. 8609 


9 


2 


00 


11 


20 


30 


41 


51 


.9320 


1.000 


. 8014 


. 8689 


3 


3 


00 


01 


10 


12 


21 


22 


1.000 


1.000 


.8571 


.9231 


4 


3 


00 


01 


10 


11 


22 


32 


.9623 


1.000 


.8400 


.8972 


-> 


o 

J 


00 


01 


11 


20 


32 


42 


» y i £- <l 


1 000 


. 8142 


. 8783 


6 


3 


00 


11 


20 


32 


42 


51 


1.000 


1.000 


.7782 


.8565 


7 


3 


00 


11 


20 


32 


42 


51 


.9722 


1.000 


.7854 


.8529 


8 


3 


00 


11 


20 


32 


42 


61 


.9510 


1.000 


.7962 


.8537 


9 


3 


00 


11 


20 


32 


42 


61 


.9361 


1.000 


. 7977 


.8496 


4 


4 


00 


01 


11 


13 


20 


32 


.9623 


.9623 


.8247 


.8747 


5 


4 


00 


01 


10 


21 


32 


43 


.9722 


.9623 


.9732 


.8497 


6 


4 


00 


11 


20 


32 


41 


53 


1.000 


.9623 


. 7928 


.8507 


7 


4 


00 


11 


23 


32 


41 


50 


.9722 


.9623 


.7955 


.8460 


8 


4 


00 


11 


20 


32 


51 


63 


.9511 


.9623 


.7901 


.8365 


9 


4 


00 


11 


20 


32 


41 


63 


.9361 


.9623 


.8009 


.8407 



5 5 

6 5 

7 5 

8 5 

9 5 


00 04 11 23 32 40 
00 11 23 32 40 54 
00 11 23 32 40 54 
00 11 23 32 40 54 
00 11 23 32 40 64 


.9722 .9722 .8078 .8560 
1.000 .9722 .7817 .8358 
.9722 .9722 .7906 .8365 
.9511 .9722 .7889 .8303 
.9361 .9722 .8034 .8385 


6 6 

7 6 

8 6 

9 6 


00 11 23 32 45 54 
00 11 23 32 45 54 
00 11 23 32 45 64 
00 11 23 32 45 54 


1.000 1.000 .7737 .8272 
.9722 1.000 .7757 .8226 
.9510 1.000 .7870 .8270 
.9320 1.000 .7930 .8278 


7 7 

8 7 

9 7 


00 11 23 32 45 54 
00 11 23 32 45 54 
00 11 23 32 45 64 


.9722 .9722 .7690 .8114 
.9511 .9722 .7847 .8203 
.9361 .9722 .7956 .8261 


8 8 

9 8 


00 11 23 32 45 54 
00 11 24 32 45 63 


.9511 .9511 .7666 .8011 
.9361 .9511 .7840 .8130 


9 9 


00 11 24 32 46 63 


.9361 .9361 .8008 .8247 
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TABLE 2E 



Table of selected GC/2 designs for block size 7 



m l 


m 2 


DESIGN 


E(10) 


E(01) 


E(ll) 


E 


2 


2 


00 


00 


01 


01 


10 


10 


11 


.9796 


.9796 


.9796 


.9796 


3 


2 


00 


00 


01 


10 


11 


20 


21 


.9796 


.9796 


.9796 


.9796 


4 


2 


00 


01 


10 


11 


20 


21 


30 


.9796 


.9796 


.9796 


.9796 


5 


2 


00 


01 


10 


11 


20 


30 


41 


.9689 


.9796 


.9280 


.9514 


6 


2 


00 


01 


10 


20 


30 


41 


51 


.9796 


.9796 


,8643 


.9236 


7 


2 


00 


10 


20 


30 


41 


51 


61 


1.000 


.9796 


.7932 


.8913 


8 


2 


00 


10 


20 


30 


41 


51 


61 


.9796 


.9796 


.7689 


.8685 


9 


2 


00 


10 


20 


30 


41 


51 


71 


.9636 


.9796 


.7963 


.8777 


3 


3 


00 


01 


02 


10 


11 


20 


22 


.9796 


.9796 


.9480 


.9635 


4 


3 


00 


01 


10 


11 


20 


22 


32 


.9796 


.9796 


.89 02 


.9286 


5 


3 


00 


01 


10 


11 


/0 


3z 


42 


.9689 


.9796 


.8585 


.9039 


6 


3 


00 


01 


10 


20 


31 


42 


52 


.9796 


.9796 


.8501 


.8990 


7 


3 


00 


10 


20 


31 


41 


52 


62 


1.000 


.9796 


.7019 


. 7956 


8 


3 


00 


10 


20 


31 


41 


52 


62 


.9796 


.9796 


.7095 


. 7953 


9 


3 


00 


10 


20 


31 


41 


52 


62 


.9636 


.9796 


. 7491 


. 8201 


4 


4 


00 


Ul 


1U 


11 




Z 3 


Jz 


.9796 


.9796 


. 8477 


.8960 


5 


4 


00 


01 


10 


11 


22 


32 


43 


.9689 


.9796 


.8214 


.8716 


6 


4 


00 


01 


10 


21 


32 


42 


52 


.9796 


.9796 


.8031 


.8568 


7 


4 


00 


11 


23 


32 


41 


52 


60 


1.000 


.9796 


.8242 


.8737 


8 


4 


00 


11 


23 


32 


41 


52 


60 


.9796 


.9796 


. 8230 


. Soil 


9 


4 


00 


11 


23 


32 


41 


52 


70 


. 9636 


. 9796 


. 8333 


. 8712 


5 


5 


00 


01 


10 


13 


20 


34 


42 


.9689 


.9689 


.8101 


.8569 


6 


5 


00 


01 


11 


23 


32 


40 


54 


.9796 


.9689 


.8273 


.8680 


7 


5 


00 


11 


21 


30 


43 


54 


62 


1.000 


.9689 


. 7785 


. 8301 


8 


5 


00 


11 


21 


30 


43 


54 


62 


.9796 


.9689 


. 8282 


. 8651 


9 


5 


00 


11 


21 


30 


43 


54 


72 


.9636 


.9689 


.8322 


.8647 


6 


6 


00 


01 


10 


24 


32 


45 


53 


.9796 


.9796 


. 7978 


.8425 


7 


6 


00 


10 


21 


33 


45 


52 


64 


1.000 


.9796 


.8087 


.8506 


8 


6 


uu 


i n 

1U 


9 1 


o o 
JJ 




jz 


04 


.9796 


.9796 


. 8268 


e% if i a 

. 8611 


9 


6 


00 


12 


21 


35 


43 


50 


74 


.9636 


.9796 


.8202 


.8524 


7 


7 


00 


12 


21 


35 


43 


56 


64 


1.000 


1.000 


.8043 


.8457 


8 


7 


00 


11 


23 


32 


45 


56 


64 


.9796 


1.000 


.8218 


.8560 


9 


7 


00 


11 


23 


32 


45 


56 


74 


.9636 


1.000 


.8232 


.8538 


8 


8 


00 


11 


23 


32 


45 


56 


64 


.9796 


.9796 


.8148 


.8464 


9 


8 


00 


11 


23 


32 


45 


56 


74 


.9636 


.9796 


.8224 


.8499 


9 


9 


00 


11 


24 


36 


45 


53 


62 


.9636 


.9636 


.8283 


.8522 
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TABLE 2F 



Table of selected GC/2 designs for block size 8 



m l 


m 2 


DESIGN 


E( 10) 


E(01) 


E(ll) 


E 


9 
4 


9 


00 


00 


01 


01 


10 


10 


11 


11 


1.000 


1.000 


1.000 


1.000 




Z 


00 


no 


m 


l n 


i ^ 

L L 


1 1 

L L 


ZU 


91 
Z L 


.9844 


1.000 


.9531 


.9746 




9 
Z 


00 


01 


10 


11 


20 


21 


30 


31 


1.000 


1.000 


1.000 


1.000 


5 


2 


00 


01 


10 


11 


20 


21 


30 


41 


.9763 


1.000 


.9606 


.9718 


6 


2 


00 


01 


10 


11 


20 


30 


41 


51 


.9746 


1.000 


.9216 


.9519 


7 


2 


00 


01 


10 


21 


30 


41 


50 


61 


.9844 


1.000 


.8877 


.9279 


8 


2 


00 


10 


20 


30 


41 


51 


61 


71 


1.000 


1.000 


.8095 


.9011 


9 


2 


00 


01 


10 


20 


30 


41 


51 


61 


.9488 


1.000 


.8771 


.9163 


3 


3 


00 


01 


02 


10 


11 


12 


20 


21 


.9844 


.9844 


.9844 


.9844 


4 


3 


00 


01 


10 


11 


20 


22 


31 


32 


1.000 


.9844 


.9163 


.9499 


5 


3 


00 


01 


10 


11 


20 


21 


32 


42 


.9763 


.9844 


.9009 


.9327 


6 


3 


00 


01 


10 


11 


20 


31 


42 


52 


.9746 


.9844 


. 8876 


.9225 


7 


3 


00 


01 


10 


20 


31 


41 


52 


62 


.9844 


.9844 


.8523 


.9006 


8 


3 


00 


11 


21 


31 


40 


52 


60 


72 


1.000 


.9844 


.8465 


.8995 


9 


3 


00 


11 


21 


31 


40 


52 


60 


72 


.9844 


.9844 


.8500 


. 8971 



4 


4 


00 


02 


11 


13 


23 


20 


32 


31 


1.000 


1.000 


.8848 


.9276 


5 


4 


00 


03 


11 


10 


23 


21 


32 


42 


.9763 


1.000 


.8751 


.9131 


6 


4 


00 


03 


11 


12 


20 


31 


42 


53 


.9746 


1.000 


.8654 


.9033 


7 


4 


00 


02 


11 


20 


31 


42 


53 


63 


.9844 


1.000 


0. 8440 


.8875 


8 


4 


00 


11 


20 


31 


42 


53 


63 


72 


1.000 


1.000 


.8473 


.8912 


9 


4 


00 


11 


20 


31 


42 


53 


63 


72 


.9844 


1.000 


.8392 


.8810 


5 


5 


00 


03 


11 


12 


21 


24 


30 


42 


.9763 


.9763 


.8695 


.9024 


6 


5 


00 


04 


10 


11 


23 


32 


41 


52 


.9746 


.9763 


.8635 


.8953 


7 


5 


00 


04 


11 


21 


30 


42 


53 


62 


.9844 


.9763 


.8617 


.8937 


8 


5 


00 


11 


23 


32 


41 


52 


64 


70 


1.000 


.9763 


.8325 


.8719 


9 


5 


00 


11 


23 


32 


41 


52 


64 


70 


.9844 


.9763 


.8271 


.8642 


6 


6 


00 


04 


11 


15 


23 


32 


40 


51 


.9746 


.9746 


.8552 


.8862 


7 


6 


00 


05 


11 


23 


32 


40 


51 


64 


.9844 


.9740 


.8450 


.8774 


8 


6 


00 


11 


23 


32 


41 


52 


64 


70 


1.000 


.9356 


.8490 


.8773 


9 


6 


00 


11 


23 


32 


41 


52 


64 


70 


.9844 


.9356 


.8295 


.8591 


7 


7 


00 


05 


11 


23 


32 


46 


50 


64 


.9844 


.9844 


.8402 


.8721 


8 


7 


00 


11 


23 


32 


46 


50 


64 


75 


1.000 


.9844 


.8155 


.8514 


9 


7 


00 


11 


23 


32 


46 


50 


64 


75 


.9844 


.9844 


.8445 


.8725 


8 


8 


00 


11 


23 


32 


46 


57 


64 


75 


1.000 


1.000 


.8461 


.8761 


9 


8 


00 


11 


23 


32 


46 


57 


64 


75 


.9844 


1.000 


.8246 


.8550 


9 


9 


00 


11 


23 


32 


46 


57 


64 


75 


.9844 


.9844 


.8369 


.8627 
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APPENDIX I 



Tables of references for cyclic designs. 

In these tables, for each initial block in the table 1, all 
other blocks of the same designs, that contain the treatment are 
given. 

Each table refers to one particular block size, ranging from 
3 to 8. Further, all blocks are classified in the same ordering as the 
table 1, that is by the number of treatments m and the reference 
letter. 
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TABLE 3A 



Table of reference for initial blocks of cyclic designs 

(block size 3) 



m 


ref 


Block 


2 




(001) (110) 


3 




(012) 


4 




(012) (013) (023) 


5 




(012) (013) (014) (023) (024) (034) 


6 


A 


(012) (015) (045) 




B 


(013) (014) (023) (025) (034) (035) 




C 


(024) 


7 


A 


(012) (014) (016) )024) (025) (034) (035) (036) (056) 




B 


(013) (015) (016) (026) (032) (045) 


8 


A 


(012) (017) (025) (035) (036) (067) 




B 


(013) (016) (023) (027) (056) (057) 




C 


(014) (015) (034) (037) (045) (047) 




D 


(024) (026) (046) 


9 


A 


(012) (015) (018) (024) (027) (035) (048) (057) (078) 




B 


(013) (014) (016) (017) (023) (025) (028) (034) (037) (038) (046) (056) 
(058) (067) (068) 




C 


(036) 
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TABLE 3B 



Table of references for initial blocks of cyclic designs 

(block size 4) 



m 


ref 


Block 


2 




(0011) 


3 




(0012) (0112) (0122) 


4 




(0123) 


5 




(0123) (0124) (0134) (0234) 


6 


A 


(0123) (0125) (0145) (0345) 




B 


(0124) (0135) (0234) (0243) 




C 


(0134) (0235) 


7 


A 


(0123) (0126 ) (0134) (0135) (0145) (0156) (0235) (0236) (0245) (0248) (0346) (0456) 




B 


(0124) (0125) (0136) (0146) (0234) (0256) (0345) (0356) 


8 


A 


(0123) (0127) (0136) (0167) (0235) (0257) (0356) (0567) 




B 


(0124) (0126) (0135) (0137) (0146) (0157) (0234) (0236) (0245) (0247) 
(0256) (0267) (0346) (0357) (0456) (0467) 




C 


(0125) (0134) (0147) (0156) (0237) (0345) (0345) (0367) (0457) 




E 


(0145) (0347) 




F 


(0246) 


9 


A 


(0123) (0128) (0145) (0156) (0178) (0246) (0247) (0257) (0348) (0357) 
(0458) (0678) 




B 


(0124) (0125) (0126) (0127) (0135) (0138) (0148) (0157) (0158) (0168) 
(0234) (0237) (0245) (0248) (0267) (0278) (0345) (0378) (0456) (0457) 
(0468) (0478) (0567) (0578) 




C 


(0134) (0137) (0146) (0167) (0235) (0238) (0256) (0268) (0347) (0358) 
(0467) (0568) 




D 


(0136) (0147) (0236) (0258) (0345) (0356) (0367) (0368) 
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TABLE 3C 



Table of references for initial blocks of cyclic designs 

(block size 5) 



m 


ref 


Block 


2 




(00011) (11100) 


3 




(00112) (00122) (01122) 


4 




(00123) (01123) (01223) (01233) 


5 




(01234) 


6 




(01234) (01235) (01245) (01345) (02345) 


7 




(01234) (01235) (01236) (01245) (01246) 
(01256) (01345) (01346) (01356) (01456) 
(02345) (02346) (02356) (02456) (03456) 


8 


A 


(01234) (01237) (01267) (01346) (01356) 
(01567) (02356) (02357) (02457) (04567) 




B 


(01235) (01236) (01247) (01257) (01367) 
(01467) (02345) (02567) (03456) (03567) 




C 


(01245) (01256) (01345) (01347) (01456) 
(01457) (02347) (02367) (03457) (03467) 




D 


(01246) (01357) (02346) (02456) (02467) 


9 


A 


(01234) (01238) (01256) (01278) (01357) 
(01456) (01458) (01678) (02357) (02457) 
(02467) (02468) (03458) (03478) (05678) 




B 


(01235) (01237) (01245) (01246) (01248) (01257) (01267) (01268) (01345) 
(01348) (01358) (01378) (01468) (01567) (01568) (01578) (02345) (02347 
(02348) (02378) (02478) (02456) (02567) (02678) (03457) (03578) (04567) 
(04568) (04578) (04678) 




C 


(01236) (01247) (01258) (01356) (01368) (01457) (01478) (02346) (02367) 
(02458) (02578) (03456) (03468) (03567) (03678) 




D 


(01346) (01347) (01367) (01467) (02356) (02358) (02368) 
(02568) (03467) (03568) 
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TABLE 3D 



Table of references for initial blocks of cyclic designs 

(block size 6) 



m 




Block 


2 




(.000111; 
















(001122) 












A 


A 


(001123) 


(001233) 


(011223) 


(012233) 








B 


Coot 223") 


Col 1233") 










5 




(001234) 


(011234) 


(012234) 


(012334) 


(012344) 




6 




(012345) 












7 




(012345) 


(012346) 


(012356) 


(012456) 


(013456) 


(023456) 


8 


A 


(012345) 


(012347) 


(012356) 


(012367) 


(012457) 


(012567) 






(013467) 


(014567) 


(023457) 


(023567) 


(024567) 


(034567) 




B 


(012346) 


(012357) 


(012467) 


(013567) 


(023456) 


(024567) 


9 


A 


(012345) 


(012348) 


(012357) 


(012378) 


(012456) 


(012468) 






(012567) 


(012678) 


(013458) 


(013578) 


(014568) 


(015678) 






(023457) 


(023478) 


(024567) 


(024678) 


(034578) 


(045678) 




B 


(012346) 


(012347) 


(012356) 


(012358) 


(012367) 


(012368) 






(012457) 


(012458) 


(012467) 


(012478) 


(012568) 


(012578) 






(013456) 


(013457) 


(013468) 


(013478) 


(013567) 


(013568) 






(013678) 


(014567) 


(014578) 


(014678) 


(023456) 


(023458) 






(023467) 


(023468) 


(023567) 


(023578) 


(023678) 


(024568) 






(024578) 


(025678) 


(034567) 


(034568) 


(034678) 


(035678) 
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TABLE 3E 

Table of references for initial blocks of cyclic designs 

(block size 7) 



m 


ref 


Block 


2 




(0000111) (0001111) 


3 




(0001122) (0011122) (0011222) 


4 




(0011223) (0011233) (0012233) (0112233) 


5 




(0011234) (0012234) (0012334) (0012344) 

(0112234) (0112334) (0112344) (0122334) (0122344) (0123344) 


6 




(0012345) (0112345) (0122345) (0123345) (0123445) (0123455) 


7 
8 




(0123456) 

(0123456) (0123457) (0123467) (0123567) (0124567) (0134567) (0234567) 


9 


A 


(0123456) (0123457) (0123458) (0123468) (0123478) (0123567) (0123578) 
(0123678) (0124567) (0124568) (0124678) (0125678) (0134568) (0134578) 
(0135678) (0145678) (0234567) (0234578) (0234678) (0245678) (0345678) 




B 


(0123467) (0123568) (0124578) (0134567) (0134678) (0234568) (0235678) 
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TABLE 3F 

Table of references for initial blocks of cyclic designs 

(block size 8) 



m 


ref 


Block 


2 




(00001111) 


3 




(00011122) (0001122) (0011122) 


4 




(00112233) 


5 




(00112234) (00112334) (00112344) (00122334) (00122344) 






(00123344) (01122334) (01122344) (01123344) (01223344) 


6 


A 


(00112345) (00123455) (01122345) (01223345) (01233445) (01234455) 




B 


(00122345) (00123445) (01123345) (01123455) (01223445) (01233455) 




C 


(00123345) (01123445) (01223455) 


7 




(00123456) (01123456) (01223456) (01233456) (01234456) (01234556) 






(01234566) 


8 




(01234567) 


9 




(01234567)(01234568) (01234578) (01234678) (01235678) (01245678) 






(01345678) (02345678) 
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APPENDIX II 



Table of reference for GC/2 designs. 

For each design in the tables 2A to 2F, in these tables all 
other blocks of the same designs that can be produced by cyclic 
development and contain treatment 00 are given. 

Each table refers to one particular block size, ranging 
from 3 to 8. Further all blocks are classified in the same ordering 
as the tables 2A to 2F, that is by the number of treatments m^ and 
m2 . 
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TABLE 4A 



Table of references for initial blocks of GC/2 designs 

(block size 3) 



TT1 

m l 


m 
m 2 


Block 


2 


2 


(00 


01 


10; 


(00 


01 


11) 


(00 


10 


11) 


3 


2 


(.00 


10 


O 1 \ 

21) 


(00 


11 


20) 


(00 


11 


21) 


4 


2 


(.00 


10 


11) 


(00 


11 


30) 


(00 


21 


31) 


5 


2 




10 


21) 


(00 


11 


40) 


(00 


31 


41) 


6 


2 


(.00 


10 


31) 


(00 


21 


50) 


(00 


31 


41) 


7 


2 


(.00 


10 


31) 


(00 


21 


60) 


(00 


41 


51) 


8 


2 


(.00 


10 


jl) 


(00 


21 


70) 


(00 


51 


61) 


9 


2 


(.00 


10 


oi \ 

3x) 


(00 


21 


80) 


(00 


61 


71) 


3 


3 




















4 


3 


(.00 


11 


11 ) 


(00 


11 


32) 


(00 


21 


32) 


5 


3 


/"aa 

(.UU 


11 


ZZ ) 


(00 


11 


42) 


(00 


31 


/ON 

42) 


6 


3 


(.00 


11 


ION 


(00 


21 


52) 


(00 


31 


42) 


7 


3 


(.00 


11 


32 ) 


(00 


21 


62) 


(00 


41 


52) 


8 


3 


(.00 


11 


32; 


(00 


21 


72) 


(00 


51 


62) 


9 


3 


1,00 


11 


32 ; 


(00 


21 


82) 


(00 


61 


72) 


4 


4 


(.00 


lz 


21) 


(00 


13 


32) 


(00 


23 


31) 


5 


4 


(.00 


12 


ft 1 \ 

21) 


(00 


13 


42) 


(00 


33 


41) 


6 


4 


(.00 


lz 


31) 


(00 


23 


52) 


(00 


33 


41) 


7 


4 


(.00 


11 


32; 


(00 


21 


63) 


(00 


42 


53) 


8 


4 


(00 


11 


32) 


(00 


21 


73) 


(00 


52 


63) 


9 


4 

*"T 


(00 


11 


32) 


(00 


21 


83) 


(00 


62 


73) 


5 


5 


(00 


12 


21) 


(00 


14 


43) 


(00 


34 


41) 


6 


5 


(00 


11 


32) 


(00 


21 


54) 


(00 


33 


44) 


7 


5 


(00 


11 


32) 


(00 


21 


64) 


(00 


43 


54) 


8 


5 


(00 


11 


32) 


(00 


21 


74) 


(00 


53 


64) 


9 


5 


(00 


11 


32) 


(00 


21 


84) 


(00 


63 


74) 


6 


6 


(00 


11 


34) 


(00 


23 


55) 


(00 


32 


43) 


7 


6 


(00 


11 


34) 


(00 


23 


65) 


(00 


42 


53) 


8 


6 


(00 


11 


34) 


(00 


23 


75) 


(00 


52 


63) 


9 


6 


(00 


14 


31) 


(00 


23 


82) 


(00 


65 


73) 


7 


7 


(00 


13 


31) 


(00 


25 


64) 


(00 


46 


52) 


8 


7 


(00 


11 


33) 


(00 


22 


76) 


(00 


54 


65) 


9 


7 


(00 


11 


33) 


(00 


22 


86) 


(00 


64 


75) 


8 


8 


(00 


13 


61) 


(00 


56 


75) 


(00 


27 


32) 


9 


8 


(00 


13 


61) 


(00 


56 


85) 


(00 


37 


42) 


9 


9 


(00 


13 


61) 


(00 


57 


86) 


(00 


38 


42) 
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TABLE 4B 

Table of references for initial blocks of GC/2 designs 

(block size 4) 



m l 


m 2 


2 


2 


3 


2 


4 


2 


5 


2 


6 


2 


7 


2 


8 


2 


9 


2 


3 


3 


4 


3 


5 


3 


6 


3 


7 


3 


8 


3 


9 


3 


4 


4 


5 


4 


6 


4 


7 


4 


8 


4 


9 


4 


5 


5 


6 


5 


7 


5 


8 


5 


9 


5 


6 


6 


7 


6 


8 


6 


9 


6 


7 


7 


8 


7 


9 


7 


8 


8 


9 


8 


9 


9 



Block 



(00 01 
(00 01 
(00 10 
(00 11 
(00 10 
(00 10 
(00 10 
(00 10 
(00 01 
(00 11 
(00 10 
(00 10 
(00 10 
(00 10 
(00 10 
(00 11 



(00 11 
(00 11 
(00 11 
(00 11 
(00 11 
(00 11 
(00 11 
(00 11 
(00 13 
(00 13 
(00 11 
(00 11 
(00 11 
(00 13 
(00 11 
(00 11 
(00 11 
(00 11 
(00 11 
(00 11 



10 11) 
10 21) 
21 31) 
21 3o) 
21 31) 
21 41) 
21 41) 
21 41) 
10 22) 

20 32) 

21 32) 
21 32) 
21 42) 
31 42) 
31 42) 
23 32) 



23 32) 
23 32) 
23 42) 
32 43) 
23 42) 
23 32) 
23 32) 
23 42) 
32 41) 
32 41) 
23 32) 
23 42) 

23 42) 
32 41) 

24 42) 
24 42) 
34 42) 
34 43) 



(00 01 
(00 11 
(00 12 
(00 11 
(00 11 
(00 11 
(00 11 
(00 02 
(00 12 
(00 11 
(00 11 
(60 00 
(00 21 
(00 21 
(00 12 



34 42) 
34 43) 



(00 12 
(00 12 
(00 12 
(00 21 
(00 12 
(00 12 
(00 12 
(00 12 
(00 24 
(00 24 
(00 12 
(00 12 
(00 12 
(00 25 
(00 13 
(00 13 
(00 23 
(00 23 
(00 23 
(00 23 



11 20) 
21 30) 
21 41) 
21 50) 
31 60) 
31 70) 

31 80) 

12 21) 

21 32) 

22 40) 
22 50) 
11 32) 

32 70) 
32 80) 
21 33) 
21 43) 
21 53) 

31 63) 

32 73) 
31 83) 
21 44) 
21 54) 
31 64) 

33 72) 

33 82) 
21 55) 
31 65) 
31 75) 

34 83) 
31 66) 
31 76) 



(00 11 
(00 10 
(00 11 
(00 10 
(00 30 
(00 20 
(00 20 
(00 12 
(00 12 
(00 11 
(00 11 
(00 21 
(00 11 



(00 11 
(00 13 
(00 13 
(00 13 
(00 23 
(00 11 
(00 23 
(00 14 
(00 14 
(00 24 



31 86) 

32 77) 

31 87) 

32 88) 



(00 14 
(00 14 
(00 15 
(00 25 
(00 25 
(00 15 
(00 25 
(00 25 
(00 15 
(00 17 
(00 16 
(00 18 



20 21) 

21 31) 

31 40) 

41 51) 

51 61) 

61 71) 
71 81) 
20 21) 

20 31) 

32 42) 

42 52) 

52 62) 
52 62) 

62 72) 

21 32) 

31 42) 

41 52) 

51 62) 

52 63) 
71 82) 

32 43) 

42 53) 

52 63) 

53 61) 
63 71) 

43 54) 
53 64) 

63 74) 

64 71) 

53 64) 
63 74) 

63 74) 

54 65) 

64 75) 

65 76) 



(00 10 
(00 11 
(00 20 
(00 31 
(00 31 
(00 41 
(00 51 
(00 11 
(00 11 
(00 21 
(00 31 
(00 31 
(00 41 
(00 51 
(00 12 
(00 22 
(00 32 
(00 32 
(00 41 
(00 52 
(00 23 
(00 23 
(00 33 
(00 44 
(00 54 
(00 34 
(00 34 
(00 44 
(00 55 
(00 35 
(00 45 
(00 55 
(00 45 
(00 56 
(00 56 



11 21) 

21 30) 
31 41) 
41 50) 
41 50) 
51 60) 
61 70) 

12 21) 

22 31) 
31 42) 
41 52) 
41 52) 

51 72) 
61 82) 

23 31) 

33 41) 
43 51) 

43 55) 

52 73) 
63 71) 

34 41) 
34 41) 

44 51) 
52 71) 
62 81) 

45 51) 

45 51) 

55 61) 
62 81) 

46 52) 

56 62) 

66 82) 
56 71) 

67 82) 
67 81) 



- 67 - 



f. 



TABLE 4C 



Table of references for initial blocks of GC/2 designs 

(block size 5) 



m 

1 


2 


Block 


2 


2 


(00 


00 


01 


10 


11) 


(00 01 


01 


10 


11) 


(00 01 


10 


10 11) 






(00 


01 


10 


11 


11) 




















3 


2 


(00 


01 


10 


11 


20) 


(00 


01 


10 


11 


21) 


(00 


01 


10 


20 21) 






(00 


01 


11 


20 


21) 


(00 


10 


11 


20 


21) 










4 


2 


(00 


01 


10 


11 


21) 


(00 


01 


10 


11 


20) 


(00 01 


11 


30 31) 






(00 


01 


10 


30 


31) 


(00 


20 


21 


30 


31) 










5 


2 


(00 


11 


20 


30 


41) 


(00 


11 


21 


30 


41) 


(00 


10 


21 


30 41) 






(00 


11 


20 


31 


40) 


(00 


11 


20 


21 


31) 










6 


2 


(00 


11 


20 


30 


51) 


(00 


11 


21 


40 


51) 


(00 


10 


31 


40 51) 






(00 


21 


30 


41 


50) 


(00 


10 


21 


30 


40) 










7 


2 


(00 


11 


20 


30 


41) 


(00 


11 


21 


30 


61) 


(00 


10 


21 


50 61) 






(00 


11 


40 


51 


60) 


(00 


31 


40 


51 


61) 










8 


2 


(00 


11 


20 


30 


51) 


(00 


11 


21 


40 


71) 


(00 


10 


31 


60 71) 






(00 


21 


50 


61 


70) 


(00 


31 


40 


51 


61) 










9 


2 


(00 


11 


20 


30 


51) 


(00 


11 


21 


40 


81) 


(00 


10 


31 


70 81) 






(00 


21 


60 


71 


80) 


(00 


41 


50 


61 


71) 










3 


3 


(00 


01 


10 


11 


22) 


(00 


02 


10 


12 


21) 


(00 01 


12 


20 21) 






(00 


02 


11 


20 


22) 


(00 


11 


12 


21 


22) 










4 


3 


(00 


01 


10 


21 


32) 


(00 


02 


12 


20 


31) 


(00 


11 


22 


30 31) 






(00 


11 


20 


22 


32) 


(00 


11 


12 


21 


32) 










5 


3 


(00 


10 


21 


31 


42) 


(00 


11 


21 


32 


40) 


(00 


13 


21 


32 42) 






(00 


11 


22 


32 


40) 


(00 


11 


21 


32 


42) 










6 


3 


(00 


11 


21 


32 


42) 


(00 


10 


21 


31 


52) 


(00 


11 


21 


42 50) 






(00 


10 


31 


42 


52) 


(00 


21 


32 


42 


50) 










7 


3 


(00 


11 


21 


32 


40) 


(00 


10 


21 


32 


62) 


(00 


11 


22 


52 60) 






(00 


11 


41 


52 


62) 


(00 


30 


41 


51 


62) 










8 


3 


(00 


11 


21 


32 


60) 


(00 


10 


21 


52 


72) 


(00 


11 


42 


62 70) 






(00 


31 


51 


62 


72) 


(00 


20 


31 


41 


52) 










9 


3 


(00 


11 


21 


32 


60) 


(00 


10 


21 


52 


82) 


(00 


11 


42 


72 80) 






(00 


31 


61 


72 


82) 


(00 


30 


41 


51 


62) 










4 


4 


(00 


02 


11 


23 


30) 


(00 


02 


13 


21 


32) 


(00 


12 


23 


31 33) 






(00 


11 


21 


23 


32) 


(00 


10 


12 


21 


33) 










5 


4 


(00 


11 


20 


32 


43) 


(00 


13 


21 


32 


43) 


(00 


12 


23 


30 41) 






(00 


11 


22 


33 


42) 


(00 


11 


22 


31 


43) 










6 


4 


(00 


11 


20 


32 


53) 


(00 


13 


21 


42 


53) 


(00 


12 


33 


40 51) 






(00 


21 


32 


43 


52) 


(00 


11 


22 


31 


43) 










7 


4 


(00 


11 


20 


32 


43) 


(00 


13 


21 


32 


63) 


(00 


12 


23 50 61) 






(00 


11 


42 


53 


62) 


(00 


31 


42 


51 


63) 










8 


4 


(00 


11 


23 


32 


50) 


(00 


12 


21 


43 


73) 


(00 


13 


31 


61 72) 






(00 


22 


52 


63 


71) 


(00 


30 


41 


53 


62) 










9 


4 


(00 


11 


23 


32 


50) 


(00 


12 


21 


43 


83) 


(00 


13 


31 


71 82) 






(00 


22 


62 


73 


81) 


(00 


40 


51 


63 


72) 
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5 
6 
7 



9 
6 
7 
8 

9 

7 



9 
8 

9 

9 



5 

5 

5 

5 

5 

6 

6 

6 

6 

7 

7 

7 

8 

8 

9 



(00 11 
(00 13 
(00 11 
(00 20 
(00 11 
(00 13 
(00 12 
(00 24 
(00 12 
(00 24 
(00 11 
(00 14 
(00 11 
(00 15 
(00 11 
(00 14 
(00 11 
(00 24 
(00 12 
(00 26 
(00 12 
(00 26 
(00 12 
(00 26 
(00 11 
(00 34 
(00 11 
(00 25 
(00 11 
(00 26 



23 34 42) 
21 32 44) 

24 32 52) 
33 44 51) 
23 34 42) 
41 52 64) 
21 34 53) 
51 63 72) 

21 34 53) 
61 73 82) 
23 34 42) 
32 43 55) 

22 34 43) 
42 53 64) 

23 54 62) 
32 43 55) 
23 34 52) 
62 73 85) 
21 34 53) 
43 55 64) 
21 34 53) 
53 65 74) 
21 34 53) 
63 75 84) 
23 36 62) 
52 63 75) 
23 35 52) 

63 74 86) 
23 35 52) 

64 75 87) 



(00 
(00 
(00 



12 23 

13 24 



13 21 
(00 13 24 
(00 12 23 
(00 33 44 
(00 14 22 
(00 42 54 
(00 14 22 
(00 42 54 
(00 12 23 
(00 24 35 
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4 4 

5 4 

6 4 

7 4 

8 4 

9 4 

5 5 

6 5 

7 5 

8 5 

9 5 

6 6 

7 6 

8 6 

9 6 

7 7 

8 7 

9 7 



(00 01 
(00 03 
(00 12 
(00 01 
(00 03 
(00 11 
(00 01 
(00 11 
(00 12 
(00 11 
(00 13 
(00 10 
(00 11 
(00 13 
(00 20 
(00 11 
(00 13 
(00 20 
(00 01 
(00 02 
(00 13 
(00 01 
(00 14 
(00 11 
(00 11 
(00 13 
(00 13 
(00 11 
(00 13 
(00 23 
(00 11 
(00 13 
(00 33 
(00 01 
(00 14 
(00 13 
(00 10 
(00 12 
(00 12 
(00 10 
(00 12 
(00 22 
(00 12 
(00 14 
(00 22 
(00 12 
(00 15 
(00 13 
(00 11 
(00 13 
(00 23 
(00 11 
(00 13 
(00 23 



10 11 22 

11 12 21 
13 22 23 

10 11 22 

11 21 32 

12 21 22 
10 21 32 
21 31 43 

13 22 
23 32 



33 
41 



20 32 42 

21 33 42 
23 32 41 
20 32 52 
31 43 52 

23 32 41 

20 42 62 

31 43 52 

10 13 20 
12 21 34 
14 21 23 

11 23 32 

22 31 42 

12 22 34 

21 30 43 

24 32 40 
24 34 43 
21 30 43 
24 32 50 
34 44 53 
21 30 43 

24 42 60 
44 54 63 
10 24 32 
21 35 42 
14 23 31 

21 33 45 

25 31 43 

22 33 45 
21 33 45 
25 31 53 

32 43 55 
21 35 43 
21 45 61 
34 43 51 
21 35 
21 36 
25 34 

23 32 

24 32 55 
34 46 55 

23 32 45 

24 42 65 
34 46 55 



43 
42 
41 
45 



23 32) 
30 33) 
30 31) 

32 43) 
30 43) 

33 43) 
42 52) 

44 53) 
40 50) 

52 60) 

53 61) 

51 62) 

52 60) 
63 71) 
61 72) 
52 70) 
73 81) 
61 72) 

34 42) 

42 43) 
33 42) 
40 54) 

43 53) 
43 51) 

54 62) 
51 61) 

51 62) 
54 62) 
61 71) 
61 72) 
54 72) 
71 81) 
71 82) 

45 53) 
43 52) 
45 52) 

52 64) 

53 64) 

51 64) 

52 64) 

60 74) 

61 74) 
50 74) 
73 82) 

65 72) 
56 64) 

54 63) 
56 62) 
56 64) 

66 71) 
61 72) 
56 74) 
76 81) 
61 72) 



(00 03 10 13 21 22 31) (00 01 12 13 22 30 31) 
(00 01 10 22 23 32 33) (00 03 13 21 22 31 32) 

(00 03 10 13 21 31 42) (00 01 12 22 33 40 41) 
(00 10 21 32 33 42 43) (00 11 22 23 32 33 40) 

(00 03 13 20 31 41 51) (00 11 22 32 42 50 51) 
(00 10 20 32 33 42 53) (00 10 22 23 32 43 50) 

(00 12 21 30 41 53 63) (00 13 22 33 41 51 62) 
(00 11 23 33 40 52 61) (00 12 22 33 41 50 63) 

(00 12 21 30 41 53 73) (00 12 33 22 41 61 72) 
(00 11 23 43 50 62 71) (00 12 32 43 51 60 73) 

(00 12 21 34 41 63 83) (00 13 22 33 51 71 82) 
(00 11 33 53 60 72 81) (00 22 42 53 61 70 83) 

(00 04 12 14 24 33 41) (00 03 10 24 32 40 41) 
(00 14 22 30 31 40 43) (00 13 21 22 31 34 41) 

(00 04 10 22 31 44 53) (00 12 21 34 43 50 54) 
(00 13 22 33 34 44 51) (00 14 20 21 31 43 52) 

(00 10 24 32 43 51 64) (00 14 22 33 41 54 60) 
(00 11 24 32 43 53 62) (00 13 21 32 42 51 64) 

(00 10 24 32 43 51 74) (00 14 22 33 41 64 70) 
(00 11 24 42 53 63 72) (00 13 31 42 52 61 74) 

(00 10 24 32 43 61 84) (00 14 22 33 51 74 80) 
(00 11 34 52 63 73 82) (00 23 41 52 62 71 84) 

(00 05 15 23 31 43 52) (00 14 22 35 43 50 51) 
(00 13 21 34 35 44 52) (00 14 21 22 31 45 53) 

(00 11 23 35 42 54 60) (00 12 24 31 43 55 65) 
(00 13 25 31 41 52 64) (00 12 24 34 45 51 63) 

(00 11 23 35 42 54 70) (00 12 24 31 43 65 75) 
(00 13 25 41 51 62 74) (00 12 34 44 55 61 73) 

(00 15 23 31 44 62 84) (00 14 22 35 53 75 81) 
(00 13 31 53 63 74 82) (00 24 40 52 61 75 83) 

(00 16 23 31 44 52 65) (00 14 22 35 43 56 68) 
(00 13 21 34 46 55 62) (00 15 21 33 42 56 64) 

(00 12 21 34 45 53 76) (00 16 22 33 41 64 75) 
(00 11 26 42 53 65 74) (00 15 31 42 54 63 76) 

(00 12 21 34 45 63 86) (00 16 22 33 51 74 85) 
(00 11 36 52 63 75 84) (00 25 41 52 64 73 86) 
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8 8 



9 8 



9 9 



(00 11 23 32 45 56 64) 
(00 13 24 32 56 67 71) 
(00 24 35 47 56 61 72) 



(00 11 23 32 45 56 
(00 13 24 42 66 77 81) 
(00 24 35 47 56 61 72) 
(00 11 24 36 45 53 62) 
(00 18 26 35 63 74 87) 
(00 37 48 52 64 73 81) 



(00 12 21 34 45 53 77) 
(00 11 27 43 54 66 75) 



74) (00 12 21 34 45 63 87) 
(00 11 37 53 64 76 85) 

(00 13 25 34 42 51 88) 
(00 17 26 54 65 78 81) 



(00 17 22 33 41 65 76) 
(00 16 32 43 55 64 77) 

(00 17 22 33 51 75 86) 
(00 26 42 53 65 74 87) 

(00 12 21 38 47 75 86) 
(00 18 46 57 61 73 82) 



- 74 - 



TABLE 4F 



Table of references for initial blocks of GC/2 designs 

(block size 8) 
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m 2 



Block 



(00 11 21 31 40 52 60 72) 
(00 10 22 31 42 51 62 70) 
(00 10 12 32 40 51 61 71) 
(00 12 20 31 41 51 60 72) 
(00 11 21 31 40 52 60 72) 
(00 10 22 31 42 51 72 80) 
(00 12 20 32 50 61 71 81) 
(00 12 30 41 51 61 70 82) 
(00 02 11 13 20 23 31 32) 
(00 02 12 13 20 21 31 33) 
(00 01 12 13 21 23 30 32) 
(00 03 10 12 21 23 31 32) 
(00 03 10 11 21 23 32 42) 
(00 03 10 12 21 31 42 43) 
(00 02 13 23 30 31 41 42) 
(00 10 21 22 32 33 41 43) 
(00 03 11 12 20 31 42 53) 
(00 01 13 20 31 42 52 53) 
(00 11 22 33 40 43 51 52) 
(00 11 21 22 30 33 42 53) 
(00 02 11 20 31 42 53 63) 
(00 13 20 31 42 52 61 62) 
(00 11 22 32 41 43 50 63) 
(00 10 21 23 32 41 51 63) 
(00 11 20 31 42 53 63 72) 
(00 11 22 33 43 52 60 71) 
(00 11 21 30 42 53 62 73) 
(00 13 21 32 41 52 63 70) 
(00 11 20 31 42 53 63 72) 
(00 11 22 33 43 52 70 81) 
(00 11 21 30 52 63 72 83) 
(00 13 31 42 51 62 73 80) 
(00 03 11 12 21 24 30 42) 
(00 01 10 13 24 31 42 44) 
(00 03 14 21 32 34 40 41) 
(00 12 20 23 31 32 41 44) 
(00 04 10 11 23 32 41 52) 
(00 01 13 22 31 42 50 52) 
(00 14 23 34 41 42 52 53) 
(00 11 23 24 30 34 42 51) 
(00 04 11 21 30 42 53 62) 
(00 10 24 31 42 51 63 64) 
(00 12 23 32 40 44 51 61) 
(00 14 21 22 33 43 52 64) 
(00 11 23 32 41 52 64 70) 
(00 14 23 34 41 52 62 73) 
(00 11 23 34 44 55 62 71) 
(00 11 21 32 44 53 62 73) 
(00 11 23 32 41 52 64 70) 
(00 14 23 34 41 52 72 83) 
(00 11 23 34 54 65 72 81) 
(00 11 31 42 54 63 72 83) 



(00 10 20 32 41 52 
(00 12 21 32 41 52 
(00 11 23 31 42 52 
(00 11 22 32 42 51 
(00 10 20 32 41 52 
(00 12 21 32 41 62 
(00 11 20 41 52 62 
(00 21 32 42 52 61 
(00 02 11 13 21 22 
(00 02 10 11 22 23 
(00 03 11 12 20 22 
(00 01 11 13 20 22 
(00 01 11 12 20 22 
(00 01 11 13 22 32 
(00 02 11 21 32 33 
(00 11 12 22 23 31 
(00 01 12 13 21 32 
(00 03 12 23 30 41 
(00 11 22 32 33 40 



(00 10 11 22 23 31 
(00 02 13 22 33 40 
(00 11 22 33 43 50 
(00 11 21 30 32 43 
(00 11 13 22 31 42 
(00 13 20 31 42 52 
(00 11 22 32 41 53 
(00 10 23 31 42 51 
(00 12 23 32 43 50 
(00 13 20 31 42 52 
(00 11 22 32 41 63 
(00 10 23 41 52 61 
(00 22 33 42 53 64 
(00 02 13 14 21 23 
(00 04 12 14 23 30 
(00 02 11 23 31 34 
(00 12 13 20 24 32 
(00 01 11 12 24 33 
(00 04 12 21 30 41 
(00 14 20 32 33 43 
(00 12 13 23 24 31 
(00 01 12 22 31 43 
(00 14 21 32 41 53 
(00 11 20 32 33 44 
(00 12 13 24 34 43 
(00 12 21 30 41 53 
(00 14 20 32 43 53 
(00 17 23 33 44 56 
(00 10 21 33 42 51 
(00 12 21 30 41 53 



61 72) 
60 70) 

62 71) 

60 71) 

61 82) 
70 80) 
72 81) 
70 81) 

30 33) 

31 33) 

31 33) 

32 32) 

33 43) 
40 43) 

40 43) 

33 44) 

43 50) 
51 52) 

41 53) 

42 53) 

51 61) 

52 61) 

52 63) 

53 64) 

61 73) 

60 73) 

62 73) 

61 71) 

61 83) 
70 83) 
72 83) 
75 85) 
32 44) 

41 43) 

42 43) 

34 43) 
42 53) 

53 54) 

44 51) 
40 54) 

54 63) 
54 60) 
54 63) 
50 61) 
64 74) 

64 72) 

65 74) 

62 74) 
64 84) 



(00 14 20 32 43 63 74 82) 
(00 12 23 43 54 6l 70 84) 
(00 20 31 43 52 61 72 84) 
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APPENDIX III 



COMPUTER PROGRAM EFFICIENCY 

This program obtains the average factorial efficiency 
for each factorial effect, and the overall efficiency. 

Input parameters to the program and their limits are: 

NF , Number of factors, less than 10 

k , block size, less than 10 

RDIM , Number of replications 

M(I) , Number of levels for each factor, less than 10 
IN(I) , the initial block. 

The limit on the size of the NN' matrix is 99x99. 
Outputs are : 

Factorial effects represented in binary form and factorial 
efficiencies . 
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JOu SATE2000 FORTRAN F1.8',60 (01) COMPILATION 1 979/0 



FORTRAN TEXT 



PROGRAM EFFICIENCY 
CCC CALCULATES THE AVERAGE FACTORIAL EFFICIENCY FO R x E A C H FACTORIAL 
COMMON NF,V,K,RDIM,NT,M( 10) ,NC(99,99),A (99, 99), AM AT (99,99) 
REAL NC 

DIMENSION S (99,99),AA(99,99),BB (99, 99 ) , W( 99 ) 
DIMENSION IN(10),MM(10) 
INTEGER RDIM,S,R,C,V 
.1000 NF=2 
K=8 
RDIM=8 
IF (NFoLE o 0) GO TO 999 
WRITE(6,4) NF,K,RDIM 

4 FORMAT (2SH PARAMETERS N,K,R,M(I) ARE : ,3 13 ) 
READ (5,20) (M(I),I=1,NF) 

WRITE (6,20) (M(I),I = 1,NF) 
20' FORMAT (1013) 
N=NF 
V=1 

DO 25 1=1, NF 
V=V*M(I) 
25 CONTINUE 
. N T = V 

WRITE(6,29) ;•>;. 

29 FORMAT (14H INITIAL BLOCK) 
READ (5,30) (I N(I) ,I=1,K) 
WRITE(6,30) (IN(I) ,I=1,K) 

30 FORMAT (1017) 
DO 40 1=1, NF 

MM (I )=M(N-I+1) 
40 CONTINUE 

DO 100 1=1, K 

LSUM=0o0 
J = NF 

L=0.0 

5 I F ( J „ E Q • 1 ) GO TO 90 

L=IN(I)/ (1 0**0-1 )) 
IN (I )=IN (I)-(L*(1 0**<J-1 )) ) 
LSUM=(LSUM+L)*MM( J-1) 
J=J-1 
GO TO 50 
90 IN(I)=LSUM+IN(I)+1 
100 CONTINUE 

WRITE (6,200) (IN(I) ,I=1,K) 
200 FORMAT (1017) 
DO 250 1=1, NT 
NC(1,I)=0.0 
250 CONTINUE 

DO 300 1=1, K 

NC(1,IN(I))=NC(1,IN(I))+1 
300 . CONTINUE 

WRITE (6,350) (NC (1, I),I=1,NT) 
350 FORMAT (99F4.0) 
CALL BLOM 
DO 399 1=1, V 
DO 397 J=I,V 



33 



JOB SATE2000 



FORTRAN Fl .B „6 Q (01) 



COMPILATION 1979/03/Uv. 



IS FORTRAN TEXT 

AM AT (I,J)=OoO 

DO 395 L=1,V 

AMAT(I,J) = NC(L,I)*NIC(L,J )+Al y lAT (.1 ,J ) 
395 CONTINUE 

AMATCI,J)=<-1)*(AMAT<I / J)/K>+1 
AMAT(J,I)=AMAT(I, J) 
IF ( I « N E .J ) GO TO 397 
AM AT <I,J)=AMAT(I,J)+RDIM 
397 CONTINUE 
399 CONTINUE 

DO 21 1=1, V 
DO 21 J = 1,V 
B(I,J) = : »0 
IFCI'.NE.J) GO TO 21 
B(I,J) =1 oO 
21 CONTINUE 
IFAIL=1 

CCC 

CCC INVERT THE V X V MATRIX AM A To PUT THE INVERSE IN Al 

CALL F04AEF (AMAT,99,B,99,V, V, A , 99 , W , A A , 99 ,BB , 99, I F A I L ) 
IF (IFAIL oNE.O GO TO 1000 

CCC 
. CCC 

CCC CALCULATE THE AVERAGE FACTORIAL EFFICIENCY FOR EACH FACTORIAL EFFE 

CALL FACTEFF 

WRITE(6,15) 
.15 FORMAT ( 1 HO) 
GO TO 1000 
999 STOP .', 
END 

CCC 
CCC 



JOU SATE20QQ 



FORTRAN F1 .!3 „6 (01) 



COMPILATION 1979/C8/C 



<ER 



FORTRAN TEXT 



SUBROUTINE BLOM 



ccc 
ccc 
ccc 



6ENERAT ES A BLOCK CYCLIC MATRIX FROM THE FIRST ROW. 

COMMON NF,V,K,RDIM,,NT,M(10)„NC(99,99),A (99, 99) , AM AT (99 99) 
REAL NC 

INTEGER RDIM,V,S,R,C 
.DIMENSION IX(10),IC(10),INC(10) 
MCNF+1) - 1 
INC(NF+1)=1 
NF1 =NF-1 
. DO 1 1 = 1, NF 
IX(I)=1 

rccr)=i 

L=N F-I + 1 

INC (L)=M(L) *INC (L+1) 
1 CONTINUE 

I R = 1 

4 I N D = N F 

5 IIX=IC(IND) 
IIY1=INC(IND+1) 
IIY2=M(1ND) 

CCC BUILDS UP A BLOCK CI RCULANT MATRIX USING IIY2 SUBMATRICES. 
CCC IR AND IIX ARE START ROW AND COLUMN POSITIONS IN THE MAIN °MATRIX 
CCC IIY1 IS DIMENSION OF EACH SUBMATRIX. ON ENTRY NC CONTAINS THE G 
CCC ROWS DEFINING THE PATTERN. ON EXIT IT CONTAINS THE BLOCK CIRCULAN 

CALL CIRC(IR,II X,IIY1,IIY2,NC) 

IC(IND) =IC( IND) +INCCIND) 

IX (IND) =IX( IND) +1 

IF ( IN D'oEQ o1 ) GO TO 100 

IFCIX(IND) .LE.M(IND-D) GO TO 4 

IX(IND) = 1 

IND=IND-1 

GO TO 5 
.100 RETURN 

END 



JOB SATE 2000 



FORTRAN F1.8.60 (01) 



COMPILATION 1 979/08/ 0'-- 



FORTRAN TEXT 

SUBROUTINE CI R C ( I R , I C , M , N, X ) 

CCC 

CCC PRODUCES A BLOCK CYCLIC PATTERN FROM FIRST GROUP OF ROWS 
CCC M IS NO„ OF ELEMENTS IN EACH ROW OF EACH CROUP* 
CCC N IS NO. OF GROUPS. 
CCC 

v COMMON NF,V,K,RDIM,NT,MD(10),NCD(99,99),A (99, 99) , AMAT (99, 99 ) 
REAL NCD 

INTEGER RDIM,S,R,C,V 
DIMENSION X(99,99) 
IRE=IR+(N-2)*M 
ICE=IC+ (N-1 )*M 
DO 100 IR1=IR,IRE,M 
IR2=IR1+M 

DO 100 IC1=IC,ICE,M 
IC2=IC1+M 

IF (IC1 iEQ .ICE) IC2=IC 
JR1 =IR1 
JR2=IR2 
JC1=IC1 
JC2=IC2 
DO 90 1=1, M 
v. DO 80 J=1,M 

X( JR2,J C2) = X( JR1,JC1) 
JC1 = JC1+1 
JC2=JC2+1 
80 CONTINUE 
JC1=JC1-M 
JC2=JC2-M 
JR1=JH1+1 
JR2=JR2+1 
90 CONTINUE 
100 CONTINUE 
RETURN 
END 



JOS S A T £ 2 1' 



FORTRAN FT .8 =6 (01) 



COMPILATION V/?v/C< ! / ! } 



ER FORTRAN TEXT 

SUBROUTINE FACTEFF 



ccc 
ccc 
ccc 
ccc 
ccc 
ccc 
ccc 
ccc 
ccc 



CALCULATES THE AVERAGE FACTORIAL EFFICIENCY, FEF, FOR. EACH EFFECT 
REPRESENTED BY AN NF-TUPLE IX(I), 1=1 /»>'. : -m ,6 AS 

FEF=(PROD*V> /DENOM , 

WHERE 

PROD = ((M (1 )**X (1) )*„. «*(M(NF)-1) **X(NF)) , 
AND DENOM=RDIM*TRACE(X*A ) , 

WHERE X IS THE CONTRAST MATRIX AND A IS A G-INVERSEo 



C 

c 

C 

c 

C 
C 



COMMON NF,V # K,RDIM,NT,M(10> ,NC( 99,9 9 ) ,A (9 9, 99 ) ,A M AT (99, 99 ) 
REAL NC 

INTEGER RDIM,S,R,C,V 
DIMENSION IXC10) 
WRITE(6,14) 

14 F0RMAT(1H0,6X,17H FACTORIAL EFFECT,2X,8H .W„EFF.) 
DO 10 I=1,NF 

:10 ix(i)=o . * 

IX(1)=1 
ALL=0.0 
20 IND=1 

IX ARRAY CONTAINS THE NF-TUPLE DEFINING THE FACTORIAL EFFECT 
CALL DF.NOM( IX , B ) 

CALCULATES DENOMINATOR OF THE EXPRESSION FOR AV FACTORIAL EFFECT-. 

PROD =1' a 
DO 30 J=1,NF 

IF (IX (J ) .NE-.O) PROD = PROD*<M<J)-1«) .. 

.30 CONTINUE 

FEF=PROD*V/B 

ALL=ALL+B " > 
WRITE(6,15) (IX(J),J=1,NF) 
i15 FORMAT(10X,1OI1 ) 

WRITE(6,16) FEF 
:16 FORMATdH + ,22X,Fl0i4) 

CHANGE IX. ORDER IS 100,010,110,001,101,011,111, FOR NF=3„ 

40 IX(IND) = IX(IND)+1 

IF( IX( I ND) „LE .1 ) GO TO 20 

IF (IND.EQ .NF) GO TO 1000 

IX ( IND) = 

IND = IND + 1 

GO TO 40 
1000 ALL=(V-1)*( V) /ALL 

WRITE (6,60) ALL 
60 FO RM AT ( 2 1 H OVERALL EFFICIENCY = ,F 6*« 4 ) 

RETURN 

END 



C 
C 

c 



I 
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FORTRAN F1 .B „6 (01) 



CORP I ;L AT 10 N 1 9 7 V / 08 / 



IER 



C 

c 
c 



FORTRAN TEXT 

SUBROUTINE DENOM(IX,B) 

COMMON NF,V,K,RDIH,NT,M(10) ,N C ( 99 ,99 ) ,A (9 9, 99 ) , AM AT (99,99) 
REAL NC 

INTEGER RDIM,S,R,C,V 

DIMENSION X (99,99 ),Y (99,99) ,IX( 10 ),Z (99,9 9) 

SET UP FIRST ROW OF FIRST MATRIX IN THE KRONECKER PRODUCT. 



C 

c 
c 



c 
c 
c 



CALL S MAT (1 ,IX,X) 

NX=M(1) 

DO 10 1=2, NF 

NY=M(I) 

1 1 = I 

CALL SMAT (II,IX,Y> 

CALCULATE KRONECKER PRODUCT o 

CALL KR ON ( X, Y , Z , NX , N Y, N Z,9 9 , 99 ,9 9) 

NX=NZ 

DO 11 11=1, NZ 

DO 11 J=1,NZ 

X(I1,J)=Z(I1, J> 
11 CONTINUE 
10 CONTINUE 

CALCULATE TRACE OF PRODUCT OF CONTRAST MATRIX AND G-INVERSE 

TRACE=0.0 
DO 15 1=1 ,V 
DO 15 J=1,V 

TR ACE=TRACE+X (I, J )*A (J, I) 
15 CONTINUE 

B=TRACE*RDIM 

RETURN 

END 



JOB SATE2C00 



FORTRAN F1 o B- o 60 (01) 



COMPILATION' 1 9 7970,*?/ G 



4ER FORTRAN TEXT 

SUBROUTINE SM AT ( 1 1 , I X, X ) 

COMMON NF,V,K,RDIM,NT,M(10) ,NC( 99,99) ,A (99,99) , AM AT (99,99) 
REAL NC 

INTEGER RDIM,S,R,C,V 
DIMENSION I X( 10),X(99,99) 
N=M (II) 

IF(IXCII) .EQ.-O) GO TO 50 
X(1,1) = M(II ) - 1 
DO 10 J=2,N 
.10 X(1,J) = -1 

CALL CI RC(1 ,1 ,1 „N,X> 
RETURN 
50 DO 60 1=1, H 
DO 60 J = 1,N 
60 X(I,J) = 1 
RETURN 
END 
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JOB SATE 2 QOO 



FOR.TR AN F1.B.60 (01) 



COMPILATION 1979/08/ 



ER 



FORTRAN TEXT 



C 

c 
c 



SUBROUTINE KR N ( X ,Y , 2 , NX ,N Y ,N Z , MX ,tt Y, KZ ) 

CALCULATES THE KRONECKER PRODUCT Z = X*Y-. 

DIMENSION X (MX,99),Y (MY,99) ,Z CMZ, 99) 
IR = 1 
IC = 1 

DO 200 1=1, NX 
DO 100 J=1,NX 
DO 90 K=1,NY 
DO 80 L=1,NY 

Z(IR,IC) = X(I,J)*YCK,L) 

IC = IC+1 
80 CONTINUE 

IC=IC-NY 

IR =IR+1 
90 CONTINUE 

IR=IR-NY 

100 CONTINUE 

I R = I R < N Y 

IC=1 
200 CONTINUE 

NZ -NX *N Y 

RETURN 

END 



